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Clusters of galaxies have not yet been detected at gamma-ray frequencies; however, the recently 
launched Fermi Gamma-ray Space Telescope, formerly known as GLAST, could provide the first 
detections in the near future. Clusters are expected to emit gamma rays as a result of (1) a 
population of high-energy cosmic rays fueled by accretion, merger shocks, active galactic nuclei and 
supernovae, and (2) particle dark matter annihilation. In this paper, we ask the question of whether 
the Fermi telescope will be able to discriminate between the two emission processes. We present 
data-driven predictions for the gamma-ray emission from cosmic rays and dark matter for a large 
X-ray flux limited sample of galaxy clusters and groups. We point out that the gamma-ray signals 
from cosmic rays and dark matter can be comparable. In particular, we find that poor clusters 
and groups are the systems predicted to have the highest dark matter to cosmic ray emission ratio 
at gamma-ray energies. Based on detailed Fermi simulations, we study observational handles that 
might enable us to distinguish the two emission mechanisms, including the gamma-ray spectra, the 
spatial distribution of the signal and the associated multi- wavelength emissions. We also propose 
optimal hardness ratios, which will help to understand the nature of the gamma-ray emission. Our 
study indicates that gamma rays from dark matter annihilation with a high particle mass can be 
distinguished from a cosmic ray spectrum even for fairly faint sources. Discriminating a cosmic ray 
spectrum from a light dark matter particle will be instead much more difficult, and will require long 
observations and/or a bright source. While the gamma-ray emission from our simulated clusters is 
extended, determining the spatial distribution with Fermi will be a challenging task requiring an 
optimal control of the backgrounds. 

PACS numbers: 98.70.Rz, 95.35+d, 98.65.-r, 98.80.Cq, 12.60.Jv 

I. INTRODUCTION 

Clusters and groups of galaxies are the largest gravitationally bound matter structures observed in the Universe. 
Although these objects are expected to host several high energy phenomena, the resulting electromagnetic non-thermal 
emission is far from being fully understood [U, • A hallmark of the occurrence of non-thermal phenomena in these 
large structures is the detection, in numerous clusters, of extended radio emission associated to the synchrotron losses 
of relativistic cosmic-ray electrons [E S S 0, 01 ■ The acceleration of cosmic rays in galaxy clusters can originate 
from a number of physical processes, including violent shocks produced in cluster-cluster mergers and the accretion 
of smaller structures B, H, M Q, the re-acceleration of cosmic rays injected by galactic sources like active galactic 
nuclei and supernovae [lO| , and inelastic collisions of primary cosmic-ray protons [l3[ producing showers of secondary 
particles, including relativistic electrons and positrons as well as gamma rays. 

The radio emission from clusters, perhaps the most solid source of observational information on non-thermal phe- 
nomena in these objects, broadly falls into two classes featuring different spatial distribution, polarization and emission 
location within the cluster. The first class, radio relics, is characterized by irregular radio morphologies and is typi- 
cally located in external regions of the cluster [3] . The second class of diffuse cluster radio sources is that of radio 
halos p^ . whose emission is typically centered on the cluster and follows a similar spatial distribution as e.g. the 
thermal X-ray emission from the intra-cluster medium (ICM) gas. Radio relics, unlike radio halos, exhibit prominent 
polarization, and are often associated to cluster regions thought to host shock activity. The origin of radio halos, 
instead, is far from being fully understood P, The upcoming generation of low-frequency radio arrays, includ- 
ing the Giant Metrewave Radio Telescope (GMRT), the LOw Frequency ARray for radio astronomy (LOFAR), the 
Murchison Widefield Array (MWA) and the Long Wavelength Array (LWA), will improve the current observational 
situation in the near future. 

In the recent past, claims of hard X-ray emission from nearby cluster of galaxies have been reported [31 . One 
possibility is that this X-ray emission originates from the inverse Compton (IC) scattering off background radiation 
of the same non-thermal high energy electron population responsible for the radio emission. Most of the clusters 
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reportedly detected at hard X-ray frequencies are merging clusters [13 , and the detections themselves are still debated 
[isl . ITgj . In addition, it is unclear if simple models for the non-thermal population fueling radio and hard X-ray 
emissions can self-consistcntly explain data from clusters (sec e.g. the case of new radio data from the Ophiuchus 
cluster of galaxies [1^ in connection with the recent claim of a non-thermal hard X-ray detection with INTEGRAL 
[2l[). Future hard X-ray missions, including the Nuclear Spectroscopic Telescope Array (NuSTAR) and the 
International X-ray Observatory (IMX) [2^, will play a key role in settling the mentioned controversies, as well as in 
fostering our understanding of processes underlying non-thermal activity in clusters . 

Gamma rays, covering the highest end of the electro-magnetic spectrum, can potentially provide essential infor- 
mation on high-energy phenomena in groups and clusters of galaxies. So far, these objects have not been detected 
in gamma-rays, and data from the EGRET telescope on board the Compton Gamma-Ray Observatory [23| have 
only produced upper limits psj (see however p6j). Statistically conclusive discoveries and a true revolution in our 
understanding of the highest energy phenomena in galaxy clusters are anticipated with the Fermi Gamma-ray Space 
Telescope, formerly known as GLAST which was successfully launched on June 11, 2008. The main instru- 
ment onboard Fermi, the Large Area Telescope (LAT), represents a remarkable jump in sensitivity compared to its 
predecessor EGRET [H]. The LAT also extends the EGRET energy range (20 MeV to 10 GeV) to much higher 
gamma-ray energies, up to about 300 GeV. Aside from probing high and ultra-high energy particle physics processes 
in astrophysical sources, the detection of gamma rays from clusters would also help to establish the properties of 
the primary proton population within clusters, and possibly clarify its role in various cluster phenomena These 
include e.g. the question of the origin of radio halos [l^, the "cooling-flow" problem [2^ and particle acceleration 
within cluster merger shocks (29j . 

Clusters of galaxies are potentially powerful observational probes of cosmology (see e.g. 113, HH HH ) . In this context, 
the accurate understanding of non-thermal phenomena in clusters is crucial to their ultimate utility as cosmological 
probes. Specifically, the estimate of cluster masses, which is at the basis of all cosmological applications using clusters, 
frequently relies on the assumption of hydrostatic equilibrium between gravitational forces and the thermal pressure 
supplied by the ICM. The accuracy of hydrostatic mass determinations is therefore limited by our understanding of 
the non-thermal pressure provided by cosmic rays, turbulence and magnetic fields in the ICM [H, [s^]. The Fermi 
Gamma- Ray Space Telescope, as well as current and planned ground based atmospheric Chcrcnkov telescopes (ACTs) 
will be able to probe the energy density supplied by cosmic rays, and possibly the evolution with redshift of cosmic 
ray pressure, to an unprecedented level of accuracy [H, [3^. 

A more exotic possibility for non-thermal activity in galaxy clusters was first envisioned by Totani in Ref. [sj: 
the pair annihilation of weakly interacting massive particles (WIMP) constituting the dark matter halo. Following 
that seminal work, Colafrancesco, Profumo and UUio calculated in |3g | the complete multi-frequency spectrum, for 
the case of the Coma cluster, resulting from dark matter annihilation. The emission spectrum extends from radio to 
gamma-ray frequencies, and includes the secondary emissions from the non-thermal electrons and positrons produced 
as final stable particles in dark matter annihilation events. In addition, Ref. [1^ also studied the heating of the ICM 
produced by the energy injected by dark matter annihilation, as well as the induced Sunyaev-Zeldovich signal. As 
far as indirect signals from particle dark matter annihilation, by far the best studied non-thermal radiative emission 
is the production of gamma-rays [s^. [40|. Specifically, when two dark matter particles annihilate, gamma rays result 
both directly from loop-suppressed diagrams as well as from the subsequent decays or radiative emission (e.g. from 
final state radiation) of standard model particles produced in the annihilation final state, like quarks, leptons and 
gauge and Higgs bosons. The resulting gamma rays have energies extending up to the kinematic limit set by the 
WIMP mass (the pair annihilation event occurs for highly non-relativistic dark matter particles), predicted to be in 



the 10-1000 GeV range in the best motivated models [4lJ, (see however [43l |44 

The LAT instrument onboard Fermi is a tremendous tool for the indirect search for particle dark matter with 
gamma rays [i^ . Of special relevance for the present study, the flux of gamma rays from WIMP dark matter 
annihilation in clusters of galaxies has been shown to be in principle large enough to be detectable by Fermi-LAT 
[ssl. [47|. Dark matter annihilation does not only produce gamma rays, but also additional stable particle species, 
such as energetic electrons and positrons. These, in turn, produce synchrotron, IC and bremsstrahlung radiation, 
with uniqu e sp ectral features. The multi-wavelength emission from dark matter annihilation was studied in detail 
in [i^, Isd . [sTI . [5^ . [ssl . [5^ . [55I . , and specifically in clusters of galaxies in [3^, . Interestingly, it was 

demonstrated that if dark matter annihilation fuels to some appreciable degree either the radio emission in the Coma 
cluster [111 , or the hard X-ray emission in the Ophiuchus cluster [131 , Fermi is almost guaranteed to have the sensitivity 
to detect gamma rays from dark matter annihilation. 

Upcoming gamma ray observations of clusters of galaxies have therefore profound implications for cosmology and, 
possibly, the discovery of New Physics. If Fermi detects gamma rays from clusters, a crucial point will be to conclusively 
assess the nature of the mechanism responsible for the emission. In the present study, we focus on how to tell apart 
gamma rays produced by standard, astrophysical mechanisms such as cosmic rays from those resulting from WIMP 
dark matter annihilation. On general grounds, we expect to have three handles to differentiate the two emission 
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mechanisms: 

• the gamma-ray spectrum: models for cosmic ray production of gamma rays predict a flux as a function of 
energy which differs from what expected out of dark matter annihilation. The low photon statistics represents 
a challenge for meaningful discrimination based on the gamma-ray spectrum. We thus study the best angular 
and energy range, and propose a technique based on hardness ratios that will help to discriminate cosmic rays 
from dark matter 

• the spatial distribution: depending upon assumptions on the dark matter substructure distribution and density 
profile, as well as on the primary cosmic ray source distribution, we predict that clusters can appear to be 
extended gamma-ray sources. We study whether this can be used to differentiate gamma rays emitted by dark 
matter from those produced by cosmic rays 

• the multi-wavelength emission: comparing the results of hydrodynamical simulations of cosmic rays in clusters 
to our predictions for a dark matter scenario, we find that the ratio of the hard X-ray to gamma-ray emission 
is a potential diagnostic to understand the origin of non-thermal phenomena in clusters. 

In our simulations and analysis, we use the latest LAT instrumental response function and observation strategy and 
the Fermi Science Tools software package which is currently being used to analyze Fermi data. 

As an application of our theoretical study, we present predictions for gamma-ray fluxes from a large X-ray limited 
sample of 130 nearby groups and clusters of galaxies. Specifically, our predictions are based on X-ray data, and on a 
fixed set of assumptions for both dark matter and cosmic rays. This also allows for a meaningful comparison of the 
dark matter to cosmic ray induced emission at gamma-ray frequencies. We present a ranking of plausible candidates 
where one might expect a bright gamma-ray signal, and of sources where the dark matter contribution is expected to 
be stronger compared to the one fueled by astrophysical cosmic rays. In particular, we discovered that low-redshift 
groups are the most promising class of objects to search for a dark matter signal from distant extra-galactic systems. 

The organization of our paper is as follows. The following Section [TT] illustrates the model we use to compute 
the gamma-ray emission resulting frtom cosmic rays pi A|) and from dark matter annihilation pi Bp . and gives 
details on the Fermi simulation setup pi C[) . Section IIIII discusses how to study the origin of gamma rays from 
clusters, including our analysis of the optimal angular region, the spectra, hardness ratios, spatial extension and 
multi- wavelength coimterparts. We present in Section llVI our predictions and ranking of nearby clusters and groups 
according to their gamma-ray emission (the complete list is provided in the Appendix). Section fVl gives a discussion 
and summary of our results, and concludes. 

II. METHODOLOGY 

In this Section we present our modeling of the gamma ray emission from cosmic rays pi Ap and from dark matter 
piBp . For definiteness, we consider the case of the Coma cluster (ssl. [sst. but what we find applies to generic low 
redshift clusters. We also discuss, in Sec. IlICi the Fermi simulation setup we employ in our study. 

A. Gamma-ray Emission From Cosmic Rays 

Several mechanisms leading to the acceleration of relativistic particles in the intra-cluster medium have been dis- 
cussed in the literature (see e.g. [1^). Most importantly, energetic arguments suggest that powerful shocks created 
in cluster-cluster mergers and in the accretion of material onto the deep cluster gravitational potential well are sig- 
nificant sources of relativistic cosmic rays H, 0, H, H, [13, [6l|, [sl, . The same shocks can also re-accelerate originally 
lower-energy particles injected into the ICM through other processes (66j . 

The common denominator to the above mentioned scenarios is that Fermi shock acceleration yields a population of 
non-thermal relativistic cosmic rays. These include primarily high-energy electrons and protons. The former efhciently 
loose energy by synchrotron emission at radio frequencies as well as through the up-scattering of background radiation 
to gamma-ray and X-ray frequencies (inverse Compton scattering). Collisions of high-energy cosmic-ray protons with 
nuclei in the ICM produce, in their hadronic debris, neutral pions promptly decaying into two gamma rays with typical 
energies and fluxes potentially observable by Fermi-LAT. In addition, the same collisions yield secondary cosmic ray 
electrons and positrons from the decays of charged pions [1, [g^, [sl] . 

At energies > 0.1 GeV and for non-merging clusters, or clusters in the intermediate or late merger stages, most 
of the gamma-ray emission is believed to stem from gamma rays produced in neutral pion decays resulting fomr the 
above mentioned inelastic cosmic-ray proton collisions Q . The energy stored in the secondary electron-positron pairs 
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is predicted to contribute at the level of 1% or less of the total power associated to the primary cosmic ray protons 
Q . In terms of the gamma-ray emission, ref . [36j estimates that only for very low average cluster magnetic fields and 
for very steep proton injection spectral indexes can the secondary inverse Compton contribution be even 10% of that 
from 7r° decay. 

Although secondary inverse Compton emission is likely subdominant compared to neutral pion decay gamma- 
ray yields, during the early stages of a merger, primary electrons can still make a significant contribution to the 
GeV radiation. The gamma-ray emission is, again, dominantly associated to the inverse Compton of, here, primary 
cosmic-ray electrons off of the microwave radiation background. The resulting inverse Compton flux is suppressed 
compared to the hadronic gamma-ray production from nuclear interactions involving non-thermal protons only as 
long as the efficiency of acceleration of hadronic species exceeds that for electrons Q. While this is indeed the 
expectation in diffusive shock acceleration theory [3|, the reader should bear in mind that the model we outline 
below might not apply to merging clusters in the early stage of a merger, or it may give underestimates of the 
cosmic ray production of gamma-rays. In general, inverse Compton emission from primary electrons is expected to 
dominate close to acceleration sites such as large scale shocks [63. Nevertheless, simulations indicate that the inverse 
Compton emission from electrons is systematically subdominant at energies relevant to the Fermi telescope compared 
to gamma-rays produced from pion decays [H, [H, [7l| . 

In this work, we consider only the typically dominant gamma-ray emission from pion decays. The inclusion of 
gamma rays from primary electrons inverse Compton scattering could affect the results of the analysis we present 
in the following ways. First, the total gamma-ray flux, when considering this additional source, will generically be 
enhanced, potentially resulting in even better prospects for the detection of galaxy clusters at gamma-ray energies 
with Fermi-LAT. Furthermore, if IC from primary electrons plays a significant role, the morphology of the gamma- 
ray emission region from cosmic rays would also be affected, leading likely to a wider extent of the emitting region 
(including for instance peripheral cluster shock regions). This could, however, potentially hinder the discrimination 
of a cosmic-ray emission from that originating from dark matter annihilation using spatial considerations (sec lIIIDj) . 
Thirdly, the spectral analysis (sec IIIIB]) will in general be affected, depending on the spectrum of the emission, and 
a broken power law feature could arise at low Fermi energies where the emission from primary electrons may become 
comparable to that from primary hadronic cosmic rays. Finally, it is possible that clusters that host a bright active 
nucleus and, concurrently, exhibit significant IC emission from an accretion shock would feature a double peak in 
the spatial distribution of high-energy cosmic-ray sources. Such a circumstance cannot be described by the radial 
power-law functional form we adopt in the simplified spatial model for the distribution of cosmic-rays described below. 

Bearing the above caveats in mind, for simplicity we assume here that the dominant source for the gamma-ray 
emission from galaxy clusters at energies relevant to Fermi originates from inelastic collisions of hadronic cosmic 
rays [s^. [6 8l . Following the arguments outlined above, we then follow the analytical cosmic-ray model outlined in 
[35l.l36l.l7l|. In this scenario, the primary proton injection spectrum is described by a simple power law, parametrized 



by a spectral index ap independent of the position in the source. In |35j a framework is outlined for incorporating the 
fireball model for very high energy cosmic ray proton interactions with the ICM as well as pion production threshold 
effects. The resulting differential source function (with units of inverse energy and volume) is 
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dE^dV, (1) 



where 

crpp = 32 (0.96 + e*-^"2.4apjj^bg^j.^ ^2) 

models the effective inelastic p-p cross section and npj is the target nucleon density in the ICM. The quantity tt-cr 
with the dimensions of the cosmic ray proton number density, has a normalization chosen so that the kinetic cosmic 
ray proton energy is proportional to the thermal energy density of the ICM. Also, ^ = 2 is the pion multiplicity, and 
dj = 0.14a~"'^'^ -I- 0.44 is a shape parameter for the 7r°-threshold physics (sH]. By integrating over all solid angles and 
dividing by 71 at and ncR^ (to be independent of a model's spatial dependence) the final differential gamma-ray source 
function (i.e. the gamma-ray flux per unit energy and unit time, per unit target and impinging cosmic-ray flux) is 
obtained. 

The spatial distribution of cosmic ray sources is modeled in terms of the ratio of the energy density in cosmic rays 
to the energy density of the thermal gas; this ratio is taken to be a power law with radius, parametrized by f3p as 



Xp{r) = XpiR,oo) { ] (3) 
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Model ap 


I3p 


CRJHC 2.1 


-0.5 


CR_HF 2.1 


1 


CR_SC 2.7 


-0.5 


CR_SF 2.7 


1 



TABLE I: Summary of parameters for cosmic ray models used. The parameter Qp indicates the primary cosmic-ray injection 
spectral index, while the coefficient f3p stands for the bias of the cosmic-ray source spatial distribution with respect to the 
cluster's thermal gas spatial distribution (see eq. [Sjl. For all models, we set the ratio of the energy density of cosmic rays 
compared to the thermal gas Xp = 0.1. 

where Xp is the ratio of the energy density of cosmic rays compared to the thermal gas, and i?5oo is the radius of 
an enclosed spherical overdensity 500 times the critical density of the universe at the source's redshift [1^1 ■ The 
parameter f3p physically reflects the possibility that the spatial distribution of the sources of high-energy cosmic rays 
deviates, via a power-law as a function of radius, from the density profile of the thermal gas in the cluster. While a 
variation in the Mach number during the course of a cluster merger can drive different injection indexes at different 
locations we neglect here, for simplicity, any spatial variation associated to Up. 

In summary, the simple model we use here depends on three orthogonal parameters: ap, that sets the spectral 
shape of both the injected primary cosmic rays and the resulting gamma rays, Xp that sets the normalization of the 
gamma-ray flux, and /3p, that (together with the gas density distribution, that can be inferred e.g. from X-ray data) 
sets the spatial distribution of the signal. We outline below a few motivated parameter space choices which define the 
benchmark models we employ to run our Fermi simulations. We then apply these models to the specific case of the 
Coma cluster. 

• As far as the injection spectral index, estimates come from theoretical arguments as well as from numerical 
simulations. For instance, Ref. Q showed that the minimum spectral index ranges from 2.1 to 2.8. In that 
range, larger values are typically associated to forward shocks and smaller masses. The larger the mass, the 
stronger the gravitational potential and the harder the predicted injection spectral index. Structure formation 
shock theory predicts injection spectral indices of 2.0 < ap < 2.5 for strong shocks [T^I, such as those expected 
for accretion shocks and strong merger shocks. We choose here ap = 2.1, 2.7 as physically motivated cases. 
Specifically, the cosmic ray spectrum in our Galaxy is observed to be a power law with ap = 2.7, which motivates 
the large ap choice. On the other hand, clusters confine cosmic rays on cosmological time scales HI] , and are 
thus expected to give rise to a harder spectrum than that of the Galaxy: ap = 2.1 is thus also a reasonable and 
motivated choice for a harder injection spectrum [s^. This choice is, in addition, consistent with the results of 
the simulations of Ref. However, if cosmic rays in galaxy clusters are accelerated in weaker merger shocks 
the expectation is one of a softer injection spectral index. In this case, strong confinement of cosmic rays in 
clusters will likely not make the spectrum much harder. Numerical simulations, though, indicate a low efficiency 
for the acceleration of cosmic rays at weak merger shocks [70| . In summary, thus, we regard our benchmark 
value ap = 2.7 as an example case rather than, strictly, as an upper limit to the cosmic ray injection spectrum. 

• The results of recent hydrodynamical simulations (see e.g. [zH) and of other recent studies [s^ motivate our 
choices of /3p = —0.5, 1 for our benchmark models. The case of (3p = —0.5, which we label as the "cuspy" profile, 
is found in simulations that include radiative effects. The profile in ([3]) with (3p = —0.5 approximates the profile 
resulting from the radiative, hydrodynamical simulations of Ref. [7l[. Non-radiative simulations indicate instead 
a flatter proflle with /3p = 1; this is an extreme scenario where the density of cosmic ray sources relative to the 
thermal gas density grows linearly with radius. Although this might be on the verge of being unrealistic, the 
effects of cooling and heating in clusters arc also somewhat uncertain, and we adopt this case as an extreme 
possibility [H, [7l| . We remark that in recent studies [t^ it has been argued that cosmic ray activity in clusters 
is dominated by turbulent rcacceleration. In this case, in the absence of shocks, the natural expectation for Pp 
would be 0. This choice falls in between the two benchmark cases we consider here. 

• Lastly, in our Fermi simulations we assume that the cosmic ray energy density is 10% of the thermal gas energy 
density, i.e. we set Xp = Xp(i?5oo) = 0.1. The choice of the normalization to the cosmic ray proton energy 
density is limited by a few constraints [s^ : 

— The gamma rays produced in inelastic collisions yielding neutral pions, which subsequently decay into two 
photons, must be consistent with the EGRET upper limits [2^ 

— The inverse Compton up-scattering of microwave photons by the secondary population produced by 
charged pion decay (as for the neutral pions these particles are produced in collisions of the primary cosmic 
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FIG. 1: Contours levels for the gamma-ray emission from cosmic rays in the Coma (left) and Virgo (right) clusters of galaxies. 
We assume that the cosmic rays have 10% of the gas energy density {Xp = 0.1), and scan the plane (/3p,Qp) where jSp 
characterizes the spatial distribution of cosmic ray sources, and ap the primary proton injection spectrum. 



ray protons with nuclei in the ICM) will give rise to non-thermal hard X-ray and soft gamma-ray emission. 
The intensity of this radiation must also be consistent with observational data. 

— Lastly, the mentioned secondary will also radiate at radio frequencies via synchrotron emission, providing 
additional constraints from radio data. 

Assuming ap ~ 2.4, Ref. shows that in the case of Coma the radio data are at the level predicted for 
Xp ~ 0.2, with appropriate assumptions on the magnetic field distribution. Tighter constraints come from high 
frequency radio data, and depend quite sensitively on ap and, more importantly, on the value of the average 
magnetic field. EGRET data, again for the case of Coma, put milder constraints {Xp < 0.45 for ap = 2.1 and 
Xp < 0.25 for ap = 2.5 (sHj). While we consider Xp — 0.1 here in our Fermi simulated observations, we will 
consider a more conservative value of Xp = 0.01 in our survey of galaxy clusters and groups in sec. IIVI and in 
the Appendix, motivated e.g. by the results of Ref. [Ill, which showed that Fermi could be sensitive (depending 
on ap) to Xp smaller than a fraction of a percent for nearby massive clusters. 

We summarize our benchmark cosmic ray models, with their names and parameters, in Table IH "CR" stands 
for cosmic rays, while "H" and "S" respectively indicate a "ifard" and a "5o/t" primary proton injection spectrum, 
ap — 2.1 and 2.7, and "C" and "F" stand for ''Cuspy" and ^^Flaf \ respectively corresponding to f3p = —0.5 and to 1. 

As pointed out in previous studies, for a wide range of parameters Fermi will be able to detect a gamma-ray signal 
produced by cosmic rays: Figure [T] shows the gamma-ray flux for the Coma and Virgo clusters for a range of ap and 
Pr,, with Xp = 0.1. In the computation, we assumed for both clusters the gas density profiles as given in Table 1 of 
[35|. The curves indicate contours of constant integrated gamma-ray fiux above 0.1 GeV, in units of 10~^ cm~^ s~^. 
For simplicity, we assume here a Fermi sensitivity for a point-source emission of 4 x 10~^ cm~^ s~^, while the EGRET 
gamma-ray limits for Coma and Virgo are given in [2^. Points below the blue line, in the central part of the panels, 
are predicted to be within the sensitivity of Fermi. The red lines indicate instead the EGRET limits: parameter space 
points below those lines are thus ruled out by current data [2^ . 

The shape of the gamma-ray emission contours is not unexpected: larger values of /3p imply smaller gamma-ray 
fiuxes, simply because they feature a flatter cosmic ray spatial distribution (see Eq. [3]), which in turn integrates to 
smaller values. In addition, we flnd that the source function in Eq. ^ implies that the largest gamma-ray fluxes 
correspond to intermediate values (2.1< ap < 2.7) of the injection spectrum. 
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B. Dark Matter Models 



The determination of the gamma-ray flux from dark matter annihilation in clusters depends on both the dark 
matter density distribution in the specific object under consideration and on assumptions on the dark matter particle 
model. For definiteness, we again consider here the case of the Coma cluster, for which a detailed study of the dark 
matter density distribution was carried out in [ssj . 

As far as the dark matter density distribution is concerned, we assume a Navarro, Frcnk and White (NEW) profile 
[T^ I for both the smooth component (with a scale radius a, to be defined in the Equations below) and for the radial 
distribution of substructure (with a biased scaling parameter a' ~ 7a, as inferred from numerical simulations of e.g. 
[tsI . [76j ^ . We adopt here the semi-analytic approach outlined in [sl] to evaluate the contribution from the smooth 
host halo and from substructures. The relevant quantity for the computation of the dark matter annihilation signal 
is a number density of particle dark matter pairs, defined as: 

In the Equation above, the first line represents the contribution from the smooth part of the dark matter halo, while 
the second line encompasses the contribution from substructures. In particular, indicates the present day mean 
matter density in the Universe, and for the function g{y), as alluded to above, we assume the NEW prescription, i.e. 

= ,2 - (5) 

xyl + x) 

The normalization parameter p' and the scale radius a can be expressed, for a given profile, as functions of the virial 
mass Mvir and of the virial concentration parameter Cvir (by "virial" we mean assuming an overdensity Avir — 343, 
see the discussion in [ll])- Following [H], we take here for the Coma cluster a virial mass Mvir — 0.9 x 10^^ M^h'^ 
and a concentration Cvir — 10. The distance to Coma is set to 95 Mpc [3^. Further, in Eq. @ above we defined a 
reference substructure normalization parameter 

(6) 



A^{a' f f -'''' dyy^ g{y) 

Finally, the substructure model is specified by the two parameters fs and A^. The first stands for the ratio of the 
total mass in subhalos over the total virial mass, 

d*^- 77r^^« = /-^^vir, (7) 

A/cut 

where dn/dM^ indicates the sub-halo mass function, and Mcut the small scale cut off mass in the matter power 
spectrum [tsI [t^ . The second term, , indicates the weighed enhancement in the number density of dark matter 
particle pairs due to subhalos. For the definition of A^ and an extensive discussion on how to assess its value we refer 
the reader to [3^, which we follow here. A^ crucially depends upon the ratio between the concentration parameter 
in subhalos over that in isolated halos at equal mass. Given a structure formation model and a dark matter density 
profile, A^ and the mentioned ratio can be traded for each other. Notice that as Eq. @ shows, A^ is not the usually 
quoted substructure boost factor. 

We employ here two sets of fg and A^ , representing a very conservative setup with a suppressed contribution from 
substructure {^^SmootK' , or "S" case) and one where instead substructures play a very significant role in setting the 
dark matter annihilation gamma-ray signal {"Boosted' , or "B" case). For the Smooth case, we assume that only 20% 
of the mass is in substructures (/s=0.2) and that the average concentration ratio of same mass host and sub-halos 
equals 2, following what quoted in [sOl- The latter assumptions yields A^ ~ 7 x 10^. In the Boosted setup, we instead 
assume /s=0.5 and a concentration ratio of 4, implying A^ ~ 7 x 10^. While smaller dark matter substructure are in 
principle possible, our present choices are realistic and compatible with the results of N-body simulations. 



^ Notice that (i) a smaller a' /a ratio is inferred from a comparison of th e g alaxy radial density profile in clusters [76l | (see also [77j'). 
and (ii) this ratio might vary from cluster to cluster (for instance, Ref. |76| finds that the ratio scales proportionally with the cluster 
concentration). We evaluated the uncertainty on the total annihilation signal for the setup we outline in this section stemming from 
considering the generous range 2 < a' /a < 10 and we find a variation of at most ±10% in the flux within one degree and of less than 
4% for the emission within 0.1 degrees, for the case of the Coma cluster. 
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Another possibility for the overaU cluster dark matter density distribution is one where the innermost profile is 
flat. This case can be physically motivated e.g. in the context of scenarios where angular momentum is effectively 
transfered between baryonic and dark matter in the process of baryon infall in the dark matter gravitational potential 
well. This process can be responsible for a significant modification to the slope of the dark matter density profile at 
small radii, leading to large core radii. In the model of Ref. [sij, the final dark matter density distribution can be 
approximated by a profile such as: 



^^-^(y) ^ (I + x)(I + x^) ' 

which we refer to as the Burkert profile [s^]. For comparison, in Fig. [10] in Sec. IIII Dl where we compare the spatial 
distribution of gamma-ray emission for different models, we also show the radial profile of gamma-ray emission from 
a Burkert profile with our two reference substructure setups, but for simplicity we do not simulate this setup. The 
radial distribution of gamma-ray flux for this profile versus an NFW profile differs only in the inner cluster regions 
{R < 0.5 degrees for Coma). 

Our choice of the particle dark matter models for the present study was motivated by considering a reasonable range 
of masses and two different dominant final state annihilation modes. In addition, the size of the pair annihilation cross 
section was fixed according to either theoretical or phenomenological arguments. In the interest of generality and in 
order to make our results easily reproducible and comparable to previous work, we do not pick specific theoretical 
particle physics frameworks, but rather we specify the dominant final state, the particle mass and its pair annihilation 
rate. This allows one to completely determine the gamma-ray emission. 

We consider a model with a relatively large mass (mwiMP = 110 GeV), a dominant W^W~ final state annihilation 
mode, and a cross section corresponding to what is expected, for that mass, for a wino-like neutralino (i.e. for the 
supersymmetric fermion corresponding to the SU(2) gauge boson), namely [av) = 1.5 x lO"^** cm^/s. We choose the 
mentioned value for the WIMP mass for two reasons: (1) we want to use the gamma-ray spectrum resulting from a 
W^W~ final state annihilation mode, which forces us to consider mwiMP > Mw ~ 80.4 GeV, and (2) we want a 
sizable gamma-ray flux, which forces us to consider a relatively light mass. The choice of 110 GeV serves both the 
purpose of avoiding fine-tuning with the W threshold and of being heavy enough to contrast it to our second WIMP 
setup choice, described below. We call this model H, for high mass, and we assume that such a WIMP has a number 
density in accord with the cold dark matter abundance thanks to either non-thermal production [83| or to a modified 
cosmological expansion at the WIMP freeze-out [s^l • 

Our second WIMP setup is a low mass model (L), featuring mwiMP = 40 GeV (the lightest mass compatible with 
grand unified gaugino masses), a dominant bb final state, and a pair annihilation cross section approximately in accord 
with what is expected for thermal production of cold dark matter, (av) = 6 x 10^^^ cm'^/s. While approximately 
n^h^ w (3 X 10~^^ cm'^/s)/((Tw) [41|, the scatter to that relation for instance in supersymmetric models (e.g. from 
resonant annihilation channels) justifies a slightly larger value, which enhances our predicted gamma-ray fluxes. 

Having specified the setup for both the dark matter density distribution and particle properties, we can compute 
the differential gamma-ray yield (number of photons per unit energy, time and surface) as the following integral over 
the line of sight: 



d$^ /■ ^, (av) ^. , ,,..dNf 



dl ^-^M^.Ur{l))-nr{E,), (9) 



where ANf/dE^{E^) stands for the differential gamma-ray yield per annihilation for final state /, as resulting from 
the Pythia [s^ Monte Carlo simulations implemented in the DarkSUSY code [86l |. 

We summarize the dark matter model parameters in Tab. [ill The resulting gamma-ray fluxes are summarized instead 
in the second column of Tab. lIIII The various gamma-ray spectra of our cosmic ray and dark matter benchmark models 
are shown in Fig. [2l 



C. Fermi Simulation Setup 



We produced simulated Fermi observations using the Fermi-LAT observation simulator tool, gtobssim, in the Fermi 
Science Tools package (v9r7) [s^l ■ We run simulations for the specific case of the Coma cluster in terms of the cluster 
distance, mass, and size. As discussed above, we choose two spectral and two spatial models for the gamma-ray 
emission from both cosmic rays and dark matter annihilation which encompass the ranges expected for these sources. 
As shown below (see Sec. IIII B|) . these models cover a wide range of gamma-ray fluxes from clusters and, therefore, 
represent a reasonable range in the possible signal-to-noise for clusters detectable by Fermi. The range and variation 
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Model 


Mass 


Final 


{av) 


/. 




ID 


mwiMP / GeV 


State 


[cm^/s] 






DMJHB 


110 


w+w- 


1.5 X 10"^" 


0.5 


7 X lO'' 


DM_HS 








0.2 


7 X 10'^ 


DM_LB 


40 


bb 


6 X 10"^" 


0.5 


7 X W 


DM_LS 








0.2 


7 X 10^ 



TABLE II: Input parameters for the dark matter models considered here. The quantity fs indicates the ratio of the total mass 
in subhalos over the total virial mass, while stands for the weighed enhancement in the number density of dark matter 
particle pairs due to subhalos. 




FIG. 2: Differential gamma-ray spectrum times gamma-ray energy squared (i.e. spectral energy density), for the various models 
considered in the present analysis. "CR" indicates cosmic ray models, while "DM" the dark matter annihilation emission. See 
the text and Tab. HI and Ull for the meaning of the other labels and the associated input parameters. 



in expected gamma-ray flux from known clusters is considered in detail in Sec. IIVI Our simulated models simply 
provide benchmarks of what could be seen for observed clusters with similar statistics. 

The simulations were run in the default scanning mode with the Pass 6 instrumental source response functions 
(P6_V1 -SOURCE). For each cluster simulation, data files defining the cluster spectrum and images defining the 
spatial distribution were fed to gtobssim (see Sec. Ill A] and III Bl for model definitions). To include the extragalactic 
diffuse background emission, we simulate an isotropic source with the power-law spectral parametrization [82] found 
in the analysis of the EGRET data [11] . We note that the Fermi background may be lower if a significant fraction of 
the extragalactic background is resolved as AGN [H, [o^, [9l|, [ol, . 

We simulate a one year observation for each combination of our benchmark spectral and spatial models for gamma- 
ray emission from cosmic rays and dark matter annihilation. We abbreviate these eight combinations as follows: 
DM_HS for a high dark matter mass and a smooth dark matter distribution, DM_LS for a low dark matter mass 
and smooth distribution, DM_HB for a high mass and a dark matter distribution with significant substructure/boost 
factor, DM_LB for a low mass and a significant substructure boost, CR_HF for a hard cosmic ray spectrum and a 
flat cosmic ray spatial distribution, CR_SF for a soft spectrum and flat distribution, CR_HC for a hard spectrum 
and a cuspy cosmic ray distribution, and finally, CR_SC for a soft spectrum and cuspy distribution. Each simulated 
source is normalized to have the total model predicted flux for the given spectral and spatial models integrated over 
a region of 10.5 degree radius (the extent of our input images to gtobssim). The input fluxes and the number of 
simulated source photons are given in Tab. IIIII Additionally we simulate five year long observations of all models to 
test the improvement to our fits with higher statistics, as discussed below. In Fig. [3] and |4] we show images of the 
five year cosmic ray and dark matter simulations, respectively, binned to one degree and including the extragalactic 
background. The range in fluxes of these models is apparent. In particular, with little boost from substructure 
(DM_HS and DM_LS) the dark matter models arc very faint. In addition, the spatial distribution of the gamma-ray 
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CR_HC CR_SC CR_HF CR_SF 



FIG. 3: Shown from left to right are images of simulations of five year Fermi observations of the Coma cluster for models 
CRJHC, CR_SC, CR_HF, and CR_SF, including the extragalactic diffuse background. Images are 20 degrees across and binned 
to have 1 degree pixels. The color bar is the same in all cases. 




DM_HB DM_LB DM_HS DM_LS 



FIG. 4: Shown from left to right are images of simulations of five year Fermi observations of the Coma cluster for models 
DM_HB, DM_LB, DMJHS, and DM_LS, including the extragalactic diffuse background. Images are 20 degrees across and 
binned to have 1 degree pixels. The color bar matches the color bar in Fig. |3] 



emission from dark matter (DM_HB) appears more extended than from cosmic rays (see also Sec. IIIIDj) . We remark 
that we neglect here additional non-galaxy-cluster point sources expected, on average, in a 20 degree field with fluxes 
comparable or larger than what we predict for our cluster emission models. Specifically, estimates of the number 
of high-latitude gamma-ray point sources based upon the extrapolation of EGRET results and/or blazar models 
[9^, or on the actual Fermi-LAT early results |96| indicate that one would expect within a 10 deg radius region 
(what we consider in our figures) ~ 1.5 sources with a flux at or above the boosted high- mass DM setup DM_HB 
(> 50 X 10-9 cm-^s-i). 

III. DEDUCING THE ORIGIN OF GAMMA RAYS FROM GALAXY CLUSTERS 

A. The Optimal Angular Region 

Prior to embarking on the spectral and spatial analysis of our gamma-ray simulations, and of addressing the 
question of the potential of Fermi to discriminate between a gamma-ray emission in galaxy clusters originating from 
dark matter versus cosmic rays, in the present Section we investigate the optimal angular cuts from a theoretical 
standpoint. We show results for our one-year simulations for all of our cosmic ray and dark matter models, as well as 
for the extragalactic diffuse cosmic-ray background. 

Given the low statistics of photon counts, we decided to include all photons with a reconstructed energy above 0.1 
GcV. Fig. [5] shows the number of photon counts inside given angular regions, specified on the x-axis, for the various 
models. The left panel shows the four cosmic ray cases, while the right panel the dark matter annihilation induced 
gamma-ray signal. In both panels, for reference, we also show our simulated extragalactic gamma-ray background. A 
naive by-eye signal-to-background inspection would indicate a small angular region (a few tenths of a degree) as the 
optimal choice. However, given the fact that we will actually be able to subtract with some efficiency the extragalactic 
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FIG. 5: The results of our simulations of one-year observations, binned by angular regions. We take as a reference case that of 
the Coma cluster, and we are summing over all photon counts in our simulations with energies above 100 MeV. The black line 
indicates the estimated extragalactic background. 



background, it makes sense to investigate the ratio of the signal to the noise, or square root of the background. In 
addition, the sheer small number of photon counts if we chose a small angular region contain very little information. 

Figure [5] shows the ratio of the signal to the square root of the background, again for all models under consideration, 
and as a function of the angular region of interest. In the most luminous cosmic ray cases (/?p = —0.5), it appears 
that the signal-to-noise is maximized for a region of interest of around one degree. For such a choice, the left panel 
of Fig. [5] tells us that we would get around 500 photons per year above 0.1 GeV inside a one degree region, which is 
approximately 50% of the overall photon flux. A choice of one degree seems thus optimal in the case of cosmic rays. 

The brightest dark matter cases, on the other hand, feature a large and "luminous" substructure content, in 
particular at large radii, making the signal significantly more diffuse than in the cosmic ray case. Given the importance 
of differentiating between a point-like and a diffuse emission from clusters, it will be important to have a large enough 
region of interest. Additionally, the signal-to- noise for the dark matter cases is maximized in regions between 1.5 and 
3.5 degrees. Again looking at Fig. [5l we confirmed that the choice of a region of interest of 3 degrees appears to be 
optimal for a dark matter type signal. 

In conclusion, we find that the analysis of the signal-to-noise and of the total photon counts would lead us to employ 
a region of interest of 3 degrees for the dark matter and of one degree for cosmic rays. In the interest of being sensitive 
to soft spectra, (the instrument Point Spread Function (PSF) would lead to the loss of most of the low energy photons 
out of an angular region of one degree) we decided to proceed with a 3 degree region for both cases. For reference, 
we mark a 3 degree radius region in figs. [3] and [H 

B. Spectral Analysis 

We examine the simulated spectra for our eight benchmark cluster models to investigate our ability to discriminate 
the gamma-ray spectra from cosmic rays versus dark matter annihilation. We extract source spectra and response 
files, including the cluster plus extragalactic diffuse background, within a 3 degree radius of the cluster center using 
the tools gtbin and gtrspgen. To account for the extragalactic diffuse emission, we extract a background spectrum 
from an annular region with an inner radius of 10 degrees and an outer radius of 12 degrees. As is clear from Fig. [51 
this outer region contains very little cluster emission. Spectra were fit using the XSPEC spectral fitting package [93| ■ 
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FIG. 6: The signal over square root of the background (S/N~ na) in simulated one-year observations, as a function of the 
angular region, for all models under consideration here. 



Model 


Flux (> 0.1 GeV) Source 
(10"^ cm-^s-i) Counts 


ftp 




mwiMP 










(GeV) 




DM_HB 


54.7 


1823 


9 r,7+U.U8 
^•^ ' -0.06 


1.51 


92.0t»6"i 


1.02 


DM_LB 


14.6 


431 


9 Q7+1-05 
^•8 '-0.49 


1.01 


< 113.5 


1.01 


DM_HS 


2.25 


73 










DMXS 


0.597 


15 










CR_HC 


37.7 


1161 


r, OQ-I-U.08 

^•'J'^-O.OT 


1.12 


24.8+1^ 


1.11 


CR_SC 


42.9 


1188 




0.95 


< 11.7 


0.97 


CR_HF 


6.76 


187 




1.01 


< 54.3 


1.01 


CR.SF 


7.71 


208 


4.05l?:i; 


1.01 


< 40.0 


1.02 



TABLE III: Summary of one year simulations and their spectral fits, for the specific case of the Coma cluster. The second 
column lists the total flux input in to the simulations (> 100 MeV), and the third column gives the total number of simulated 
cluster photons (all angles). Columns 4 and 5 give the best-fit slope and reduced for the CLUSTERCR fits, while columns 
6 and 7 list the best-fit particle mass and reduced for Ihe DMFIT fits. Upper/lower limits on spectral parameters refer to 
90% confidence limits; all other errors are 1 sigma. Reduced chi-squared (Xv) in bold indicate a fit probability of less than 1% 
and in italics indicate that the best fit is found for the DMFIT lower mass limit of 10 GeV. 



Within XSPEC, the background spectrum is subtracted from the source spectrum after adjusting the BACKSCAL 
header keyword in the background file to account for the difference in area between the background spectrum and 
the source spectrum. 

XSPEC allows one to include custom user models, and we use this feature to include models for the gamma ray 
spectrum from both cosmic rays and dark matter annihilation. For the dark matter spectrum, wc use the routine 
DMFIT, which wc presented in ref. [1^. DMFIT is a tool that provides the gamma-ray flux from generic WIMP 
pair annihilation (i.e. from dark matter particles with specified mass and branching ratios into Standard Model final 
state annihilation modes). DMFIT is based on the same set of Monte Carlo simulations used in DarkSUSY [s^ and 
incorporates a wide variety of annihilation modes. Two data files contain the Monte Carlo simulation results giving 
the differential and integrated gamma-ray fluxes at given energies. The simulation results are then interpolated given 
the dark matter particle mass and annihilation final states supplied by the fitting routine. We additionally include 
the e~^e~ channel, where gamma-rays are radiated in the final state via internal bremsstrahlung. The e~^e~ channel - 
presently not included in the DarkSUSY code - is relevant for various non-supersymmetric WIMP models, including 
the Kaluza-Klein dark matter of Universal Extra Dimensions 42| (see Ref. [93| for model-independent limits on the 
annihilation cross section of dark matter to e^e~). 
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FIG. 7: Shown are fits to the one year simulated spectrum for model DMJHB (mwiMP = 110 GeV), for the specific case of 
the Coma cluster. The left panel shows the best-fit dark matter model using DMFIT with a W^W" final state. The right 
panel shows to same simulated data and best-fit cosmic ray spectrum using CLUSTERCR. The cosmic ray model is ruled out 
at better than 99.9% confidence. 



DMFIT consists of two data files and one Fortran routine, and the code is publicly available fi-om the authors upon 
request. DMFIT essentially reverse-engineers the use of the DarkSUSY package for the computation of gamma-ray 
spectra: given an observed gamma-ray spectrum, DMFIT allows one to fit for the best matching particle dark matter 
mass, its pair-annihilation rate and its branching ratios. In conjunction with virtually any fitting package, like XSPEC 
and the gtlike routine in the Fermi Science Tools, DMFIT can be used to reconstruct confidence level ranges for the 
mentioned particle dark matter properties. 

While the Pythia Monte Carlo simulations extend down to a WIMP mass of 10 GeV, DMFIT allows one to 
extrapolate to l ower masses. Very light WIMPs have been recently shown to be relevant even in the conte xt of 
supersymmetry jlOOl |. and they can possibly play a role in explaining the puzzling DAMA/LIBRA signal |lOl| . We 
note that in the current XSPEC version of DMFIT, the default lower limit on the dark matter particle mass is 10 
GeV. For the spectral fits shown below with only upper limits on the particle mass, the best fit typically saturates at 
the 10 GeV limit. 

The spectrum from cosmic ray interactions is modeled with the simple analytic form given in Eq. ([1]) , and included 
in XSPEC using a Fortran routine as a model named CLUSTERCR. Both XSPEC models are available from the 
authors upon request, and both models are currently being incorporated into the Fermi-LAT Science Tools as part of 
the gtlike likelihood fitting tool. 

Tab. HIT] summarizes the spectral fits to simulations of one year observations of the Coma cluster. For each simulated 
cluster spectrum, we fit both a cosmic ray and a dark matter model, regardless of the input model, to compare how 
well we can distinguish these two scenarios based on the spectrum. For simplicity, in the case of the cosmic ray 
simulations, we assume a bb final state for the DMFIT spectral model, while for the dark matter simulations, we use 
the dominant final state {bb for DM_ L and iy+ iy~ for DM_II). The overall spectral shape of the bb final state is quite 
similar to most other final states [86l . ll02l . Il03| . including W^W^ , but we refer the reader to ref. [i^ for a discussion 
of the systematic affects of the assumption of final state on the reconstruction of dark matter particle properties. 
Columns 4 and 5 of Tab. IIIII give the best-fit slope and reduced for the CLUSTERCR fits, while columns 6 and 7 
list the best-fit particle mass and reduced for the DMFIT fits. 

For a higher dark matter mass and a reasonably bright source, as is the case for DM.HB (mwiMP = HO GeV), we 
find a good fit to a dark matter spectrum, but we cannot get a good fit to a cosmic ray spectrum. Fig. [7] shows a 
comparison of the best-fit DMFIT (left panel) and CLUSTERCR (right panel) models for this simulation compared 
to the data. In this case, the spectrum is only consistent with a dark matter interpretation; the best-fit dark matter 
particle mass is low by ~ 15% but consistent with the true mass within a couple of sigma. On the other hand, our 
low dark matter mass model DM_LB (mwiMP = 40 GeV), which is also significantly fainter, can be well fit by either 
a cosmic ray spectrum with a fairly steep slope or with a dark matter spectrum. 

A similar result is seen when considering the the simulated cosmic ray models. In general, these are consistent with 
either a cosmic ray spectrum or a dark matter spectrum with a low mass dark matter particle. The cuspy cosmic ray 
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Model 


Flux 


(> 0.1 GeV) Source 
^ cm~^s~^) Counts 






n^WIMP 






(10- 






(GeV) 




DM_HB 




54.7 


9171 


2.10+"o-o2 


2.94 


100.5+:^-^ 


0.96 


DM_LB 




14.6 


2343 




1.03 


31.7ti« 


0.97 


DM_HS 




2.25 


368 


> 2.19 


0.98 


< 200 


1.00 


DMXS 




0.597 


100 










CR_HC 




37.7 


5819 


9 90+U.U3 


1.11 


29.2+f^ 


2.02 


CR^C 




42.9 


5990 


9 7C+0.06 
^- ' "-0.06 


0.99 


< 11.1 


1.40 


CR.HF 




6.76 


977 


2.66l°J? 


0.99 


24 1+**° 


1.03 


CR_SF 




7.71 


1145 


o 09+0.25 


0.92 


< 13.0 


0.98 



TABLE IV: Summary of five year simulations and their spectral fits. The second column lists the total flux input in to the 
simulations (> 100 MeV), and the third column gives the total number of simulated cluster photons (all angles). Columns 4 
and 5 give the best-fit slope and reduced for the CLUSTERCR fits, while columns 6 and 7 list the best-fit particle mass and 
reduced for the DMFIT fits. Upper/lower limits on spectral parameters refer to 90% confidence limits; all other errors are 
1 sigma. Reduced chi-squares in bold indicate a fit probability of less than 1% and in italics indicate that the best fit is found 
for the DMFIT lower mass limit of 10 GeV. 

distribution gives fluxes 5-10 times higher than a flatter distribution, and in these cases the possible dark matter mass 
is limited to be quite low (24.8l4'g GeV for a hard cosmic ray spectrum and a 90% upper hmit of 11.7 GeV for a soft 
cosmic ray spectrum). Even for the fainter CR_HF and CR_SF models, we find a 90% upper limit on the possible dark 
matter mass of 40 — 55 GeV. Interestingly, the cosmic ray spectral slope, ap, is overestimated in all cases, though the 
best fits are within a couple sigma of the true slope. This is most likely due to the hard spectrum of the extragalactic 
background emission. As the dark matter particle mass is also somewhat underestimated for DM_HB, it appears that 
the background subtraction leads to a general softening of the source spectrum. The effect of the background can be 
mitigated and the cosmic ray spectral slope recovered to good accuracy if a smaller energy rage extending only up to 
10 GeV is used in the spectral fits. However, with such a low high energy cutoff, a high mass dark matter particle 
model, like DM_HB, is indistinguishable from a hard cosmic ray spectrum. We, therefore, recommend using a large 
energy range {r^ 0.1 — 150 GeV) to investigate the dominant source of the emission but a lower high energy cutoff 
(~ 0.1 — 10 GcV) to determine the model parameters. 

As Fermi is expected to have at least a five year mission lifetime, we also consider the improvement in the spectral 
constraints for five year simulations of the same models; the results are shown in Tab. IIVI With deeper observations, 
even a low mass dark matter particle model is ruled out for the brighter cosmic ray models, CR_HC and CR_SC, and 
for the fainter models, CR_HF and CR_SF, the upper limit on the particle mass is significantly decreased requiring a 
very light dark matter particle. As an example, we show in Fig. [8] the best-fit DMFIT (left panel) and CLUSTERCR 
(right panel) models to the five year simulated spectrum of CR_HC {ap = 2.1). A dark matter model is clearly 
inconsistent with the spectral data. As noted above, here we have only considered a bb final state, a reasonable 
assumption as it has a similar spectral shape to most other potential final states. However, we note that if we fit a 
T+r~ final state, which has a significantly harder gamma-ray spectrum than bb, we get a much worse fit to the cosmic 
ray simulations. With five years of data, the best-fit cosmic ray slopes for the CLUSTERCR model are closer to the 
input values, but they are still too high by ~ la for the soft spectral model and ~ 3a for the hard spectral model. 
Again, this offset can be removed by using a smaller energy range (0.1-10 GcV) in the spectral fit to reduce the effects 
of the background. 

For the dark matter simulations DM_HB and DM_LB, the longer exposure time significantly reduces the errors 
on the particle mass estimates, decreasing both the systematic shift to lower masses and the statistical errors. The 
simulated high particle mass spectrum is again very inconsistent with a cosmic ray-type spectrum. However, the low 
mass model with its moderate fiux is still well fit by a cosmic ray model. 

Unfortunately, the simulations of emission from dark matter annihilation with a smooth dark matter spatial distri- 
bution (i.e. little substructure) have very low fluxes and low S/N even with five years of data. For these simulations 
(DM_HS and DM_LS), we cannot get good constraints on the spectral model. The best case with the smooth spatial 
distribution is the five year simulation of DM_HS, and as can be seen in Tab. llVi we can only derive a lower limit on 
the possible cosmic ray slope and an upper limit on the possible dark matter particle mass. 

In summary, assuming a dominant cosmic-ray or dark-matter origin, gamma ray emission from dark matter anni- 
hilation with a high dark matter particle mass (mvviMP > 50 GeV) can be distinguished from a cosmic ray spectrum 
even for fairly faint sources (one year data for CR_HF and CR_SF, for example). Distinguishing a cosmic ray spectrum 
from a low dark matter particle mass is much more difficult and requires deep data and/or a bright source. 

While our models show that either cosmic rays or dark matter annihilation can dominate the gamma-ray emission 
from clusters, a mix of the two is another likely scenario. As a final test, we simulate two clusters with emission from 
both cosmic rays and dark matter with equal fluxes for five years of observing time. The first simulation combines 
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FIG. 8: Fits to the five year simulated spectrum for model CR_HC (ofp = 2.1). The left panel shows the best-fit dark matter 
model using DMFIT with a bb final state, and the right panel shows the best-fit cosmic ray spectrum using CLUSTERCR. 
The dark matter model is ruled out at better than 99.9% confidence. Note that in the left panel data points above the best-fit 
particle mass are not shown as the fiux here is predicted to be zero. 



emission from models DM_HB and CR_SC each normalized to have a flux of 2.5 x 10"'* photons m~^ s~^; the second 
simulation combines emission from DM_LB and CR_HC each with a fiux of 1.3 x 10"''^ photons m"^ s"*. Note that 
the total cluster fiuxes from both sources of emission are then within the range of the original models. The first case 
is more optimistic combining a high mass dark matter particle with a soft cosmic ray spectrum, while the second case 
has about half the total fiux and combines a low particle mass with a hard cosmic ray spectrum. We do not simulate 
fainter mixed models, as we have already demonstrated the difficulty of distinguishing the type of emission for faint 
sources. 

For the more optimistic case of a high mass dark matter particle and a soft cosmic ray spectrum, we find that 
neither a cosmic ray only spectral model (xi/ = 1.35 for 94 DOF) nor a dark matter only model {xv — 2.03 for 94 
DOF) provide a good fit. Fitting to a combined dark matter and cosmic ray spectral model, however, gives a good fit 
(Xf = 0.87 for 92 DOF) with the normalizations of both components non-zero at > 3cr. Here the DM particle mass 
is underestimated and the cosmic ray slope overestimated by 2 — 3cr with mwiMP = 87l^ GeV and ap = 3.5 ± 0.4. 
If the dark matter particle mass is known from other astronomical observations or direct detection experiments, then 
the errors can be significantly reduced. Fixing the dark matter mass at the true value, we find ap = 3.1^g'2 and the 
normalizations of the cosmic ray and dark matter components are both determined to ~ 15%. The results for the 
second mixed simulation with a low particle mass and hard cosmic ray spectrum unfortunately show that one cannot 
distinguish the source of the emission. Comic ray only, dark matter only, and mixed models all give good fits to the 
data even with the dark matter particle mass fixed at its true value. 

Overall, we find, for bright enough sources, that if the dark matter contribution to the cluster flux is significant and 
the particle mass is not very low (rnwjMP > 40 GeV) the presence of a dark matter component can be seen even in 
the presence of a significant gamma-ray flux from cosmic rays. The reconstruction of the model parameters is more 
difficult in this case, but can be significantly improved if something about the dark matter mass in known. 

C. Hardness Ratios 

In addition to the spectral analysis, we also investigated the use of simple energy band ratios, hardness ratios, 
to differentiate the gamma-ray signal of cosmic rays and dark matter. Hardness ratios give a quick estimate of the 
slope/shape of the spectrum and can be used as a rough spectral analysis for lower photon count data (for instance, 
dimmer sources or shorter observing time) when it is not possible to accurately estimate the flux in more than a 
couple of energy bins. 

We defined two energy bands hard (H) and soft (S), and our hardness ratio is defined as HR = (H - S)/(H + S). 
A couple of considerations were used in defining these bands. We chose a lower energy limit of ^ 150 MeV so that 
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FIG. 9: The hardness ratio for cosmic ray and dark matter models. Indicating with H the total number of (hard) photons 
with energies above 0.7 GeV, and with S those (soft) in the energy band between 0.15 and 0.7 GeV, we define here as hardness 
ratio the quantity (H-S)/(H+S). The top plots are for high flux point sources without background, used as reference results. 
The middle and bottom rows of plots show the results for the simulated Fermi cluster observations of Coma, including the 
uncertainties due to the background subtraction. The red line in each plot is the same as in the top plots (just connected, not 
a best fit). 



the expected cosmic ray spectrum is a simple power law (e.g. away from threshold effects). We also require that the 
hard band have at most a factor of 10 difference in photon counts from the soft band for the range of models under 
consideration. The energy bands we chose are S = 0.15 < i? < 0.7 GeV and H = > 0.7 GeV. This particular choice 
of bands has the additional convenient feature that the expected ratio, HR, is negative for cosmic rays and positive 
for dark matter. 

As for the spectral analysis of our Fermi simulations, we used a 3 degree radius source region and an annulus from 
10 to 12 degrees to estimate the background. The photon counts in the background region are subtracted from the 
counts in the same energy band in the 3 degree source region, after accounting for the difference in area. We compared 
this annular background signal to a simulation of the background in the source region, and they match very well. 

To find the expected hardness ratios for the different spectra after including the Fermi instrument response, we 
ran simulations of point sources for cosmic rays {ap between 2 and 3) and dark matter (neutralino mass between 40 
and 400 GeV, b-b final state), without any background and set to very high fluxes (top panel of Fig. For all of 
these simulations, we took 10^ photons. Here we indeed see that the expected ratio, HR, is negative for cosmic rays 
and positive for dark matter. Also, at least with such good statistics, the ratio correlates well with ap or the dark 
matter mass. The results of our cluster simulations agree well with these results. The hardness ratios for our eight 
simmulated cluster models, after background subtraction as described above, are plotted in the middle and lower 
panels of Fig. [9l for 1 and 5 years of data. 

All of these simulations agree well within errors with the results from the simple high-flux analysis, though the errors 
are large for the fainter cases. While it can be difficult to tell, for example, a lower dark matter mass from a harder 
cosmic ray spectrum, these results show that it is possible to use this simple two energy band ratio for distinguishing 
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FIG. 10: The gamma-ray emission as a function of the angular distance from the center of the Coma cluster, for cosmic rays 
(left) and for dark matter (right). The first two lines in each panel correspond to our reference setups, and in both panels the 
thin solid black lines offer a comparison for dark matter (left) and cosmic rays (right). 



some cosmic ray and dark matter models. The hardness ratio also correlates well with spectral index Up or the dark 
matter mass. So, for instance, if it is known that a gamma-ray source is mostly dark matter in origin, this ratio 
provides an indication of the particle mass. From comparing the simulations both with and without background we 
find that significant errors come from the uncertainty in the background and its subtraction. Accurate modeling of 
the diffuse gamma-ray background could significantly improve this analysis. 

D. Spatial Extension of the Gamma-ray Emission 

In addition to their spectra, the gamma-ray emissions from cosmic rays and dark matter annihilation in clusters 
are expected to have different spatial distributions. With sufficient signal, cluster emission from cither source may 
also be detectable as extended, which could distinguish a cluster signal from the emission from e.g. a nearby AGN. In 
Fig. 1101 we show the gamma-ray flux versus angular distance from the cluster center for a range of models including 
our benchmark simulated models. The left panel of Fig. [TO] shows the radial distribution of fiux for cosmic ray models 
with varying Pp. Obviously, as Pp increases the fiux falls off less steeply with radius, but gcnerically all of the cosmic 
ray models are relatively flat within 0.1 degrees and then fall-off steeply with increasing radius. In comparison, the 
emission from the dark matter models (right panel of Fig. [TO]) has a much flatter radial distribution, particularly 
when a significant substructure component is present. With a smaller contribution from substructure, as is the case 
for DM_IIS, the gamma-ray flux falls steeply at small radii but eventually flattens around R ~ 0.3 degrees where 
the contribution from substructure becomes signiflcant. For comparison, we also show in the right panel of Fig. [TO] 
an alternative dark matter distribution to the centrally cuspy NFW model, the Burkert profile discussed in Sec. Ill Bl 
(labeled Burk. in the figure). The radial profile of emission from the Burkert density distribution is very similar to 
our simulated NFW model at large radii, but flatter at small radii {R < 0.5 degrees). In practice, this difference can 
only be resolved for bright sources and at high energies where the PSF is smaller. 

We now turn to our simulations and ask whether our simulated models can be detected as extended and whether 
a difference in spatial distribution can be observed for our cosmic ray versus dark matter models. Unfortunately, a 
full joint fitting of the spectral and spatial distribution of a source cannot currently be carried out with the Fermi 
Science Tools, and here we choose to separate the spatial an spectral modeling. For simplicity, we consider only 
one simulated spectrum (plus the two reference spatial models) each for the dark matter and cosmic ray cases, the 
high mass spectrum for dark matter and the soft spectrum for cosmic rays. We also fit only the higher S/N, 5 year 
simulations. For each of the four simulations considered (DM_HB, DM_HS, CR_SC, CR_SF), we create images of the 
cluster plus extragalactic background for E > 100 MeV and R < 3 degrees binned to have 0.2 degree pixels. This 
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TABLE V: Results of spatial fitting to the 5 year simulations. Columns 3, 4, and 6 list the value of the Cash statistic for 
the fit over the number of degrees of freedom. Values of the Cash statistic in bold indicate fits with a probability of less than 
1%. Column 2 indicates whether or not the extragalactic diffuse background was included. In row 2, for the /3-model fit to 
DM_HB including the gamma-ray background, a pegs at its upper limit and the upper limits on the parameters of the model 
are unconstrained. 

binning was chosen to give a reasonable number of counts per bin while still being smaller than the Fermi-LAT PSF 
at all but the highest energies. While the PSF can be significantly reduced by only considering higher energies, we 
did not wish to reduce the photon counts from already faint sources. 

First, we create a model of the energy averaged PSF for each spectral model considered by simulating a high flux 
point source with that spectrum. An image of the point source was then created with the same binning, energy cut, 
and angular region as for our cluster simulations. This image, re-normalized to one, is used as a PSF model, which 
is convolved with a given spatial model and then fit to the simulated cluster data. All spa tial fitting is done using 
the package Sherpa |l04l |. distributed as part of the Chandra data analysis software CI AO [lOSj . We consider three 
different models for the spatial distribution as provided by Sherpa: a delta function to test if a source is point-like, a 
Gaussian, and a /3-model of the form 

S{r)^So(l+(^y] . (10) 

Neither a Gaussian nor a /3-model exactly describes the spatial distributions shown in Fig. (TUl but these models give 
an indication of t he ex tent and slope of the source distribution. We fit using the maximum-likelihood based Cash 
statistic, CSTAT |l06l |. which is more appropriate than a statistic for data with few counts per bin but has the 
property that AC is distributed approximately as Ax^ when the number of counts in each bin is > 5 (as is the case 
for our images). Our results are summarized in Tab. |Vl 

Initially, we consider only the simulated cluster emission and neglect the extragalactic background. Tab. |V] reveals 
a couple of trends. First, none of the simulated clusters are consistent with a point source, with the exception of the 
very faint model DM_HS whose spatial distribution is not well constrained. Second, as expected, the dark matter 
models are more extended, in terms of the Gaussian FWHM, and have flatter profiles, in terms of the /?-model slope 
a than the cosmic ray models. Typically both the Gaussian and the [3 models give acceptable fits to the data. 

Finally, we consider how well we can model the spatial distribution if the extragalactic background is included. Here 
again we create images of our simulations with the same radius, binning, and energy range, but with the extragalactic 
diffuse emission included. As in the spectral analysis, we use an outer annulus between 10 and 12 degrees to measure 
the background level. The average background per pixel from this region is then added as a constant to our spatial 
modeling. For the two fainter simulated clusters, DM_HS and CR_SF, unfortunately, the addition of the background 
means that the spatial model is not well constrained; the fits are poor and the errors on the model parameters are 
large. For the brighter simulations, DM_HB and CR_SC, we again find a much better fit to an extended source than 
to a point source. These fits are noted in Tab. |Vl The fits typically worsen somewhat in terms of fit probability 
given the imperfect modeling of the background, but in general, the fit parameters such as the FWHM and slope are 
consistent within the errors with the no background results. 

In summary, clusters are expected to be extended gamma-ray sources, and with the inclusion of the contribution 
from substructure to the dark matter distribution, the gamma-ray emission from dark matter annihilation is predicted 
to be flatter and more extended than for cosmic rays. However, these differences can only be detected, given the 
extragalactic background, for bright sources with deep data. If the extragalactic background used in this work is 
an overestimate of the Fermi background, due for example to AGN unresolved by EGRET, then the situation could 
improve significantly. On the other hand, with real observations we will also have the additional uncertainty of 
imperfect knowledge of the spectrum. As a final note, what we have presented here is a relatively simple approach to 
the spatial modeling of Fermi data based on tools currently available. This analysis could be significantly improved 
through the development of tools to jointly fit the spectral and spatial distribution to arbitrary functions. 
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FIG. 11: The ratio of the gamma-ray to hard X-ray integrated flux for dark matter, as a function of the dark matter particle 
mass, for average magnetic fields B — 3 nG (left) and B = 1 fiG (right). The red arrows indicate the range predicted for 
cosmic rays as assessed from the numerical simulations of [lOTi ], with various assumptions on the cosmic ray physics. 



E. Multi-wavelength counterparts 



Several recent studies highlighted the importance of secondary radiation emitted from electrons and positrons 
produced in WIMP pair annihilation as a powerful indirect dark matter diagnostic. The non-thermal population of 
light stable leptons produced in dark matter halos via particle annihilation radiates, in fact, through synchrotron, 
inverse Compton and bremsstrahlung off the ICM gas, with peculiar spectral features. First discussed in [s^l for the 
case of galaxy clusters, the multi- wavelength emission from dark matter annihilation was also studied in detail in [51 1 
for galactic dark matter clumps and in [sl] for the dwarf spheroidal galaxy Draco. Ref. [s^ extended those analyses 
and studied constraints on particle dark matter properties from X-ray observations of nearby dwarf galaxies. Other 
recent studies include an interpretation of the significant non-thermal X-ray activity observed in the Ophiuchus cluster 
in terms of IC scattering of dark-matter produced e~^e~ [IJ, an analysis of the broad-band dark-matter annihilation 
spectrum expected from the Bullet cluster [57| and from the galactic center region [sij . In addition, radio emission 
from e+e~ produced in dark matter annihilation was considered as a possible source for the "WMAP haze" in the 



seminal paper of [55| , and subsequently analyzed in detail in [56[ , [108j and |109j . Other studies [48|, |49|, l50[ have also 
addressed synchrotron radiation induced by dark matter annihilation. 

It was pointed out in [i^] that a unique aspect of the multi-wavelength spectrum from WIMP dark matter annihila- 
tion in clusters is a strong emission at hard X-ray frequencies {E^ > 10 keV), with strikingly uniform spectral features, 
independent of the specific WIMP model. Other studies, in particular [5^, showed that even in systems where cosmic 
ray diffusion and leakage play an important role, like dwarf galaxies, X-ray and gamma-ray constraints on particle 
dark matter with current telescopes are comparable. It is therefore important to assess whether information can be 
extracted by comparing the X-ray and the ga mma- ray emission from dark matter with the results of hydrodynamical 
simulations of cosmic ray physics in clusters [107| |. The latter hard X-ray emission is thought to result from the 
inverse Compton up-scattering of mostly cosmic microwave background photons by non-thermal cosmic ray electrons 
accelerated to relativistic energies in the ICM through the multiple mechanisms summarized in the Introduction. 

While a dedicated analysis of the spectral features expected in the hard X-ray band for cosmic rays and for dark 
matter is beyond the scope of the present analysis, we wish here to compare the cumulative flux in these two frequency 
regions. Specifically, we consider the ratio of the integrated gamma-ray flu x (for E > 0.1 GeV) to the integrated hard 
X-ray flux (for i? > 10 keV) as resulting from the numerical simulations of [l07j |. and quoted in Table 3 there. Fig. [TT] 
indicates the results of [107| with red arrows. In the left panel we employ a scaling m agne tic field value of 3 /xG, whil e 



107| . As specified in ref. [10; 



in the right panel of 10 fiG, with the spatial distribution assumptions outlined in 
where we refer the Reader for further details, simulation S3 includes thermal shock heating, radiative cooling, star 
formation. Coulomb and hadronic cosmic ray losses, and cosmic rays from shocks and supernovae. Simulation S2 
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neglects the latter component, while SI also neglects radiative cooling and star formation. 

We also show in Fig. [11] the results of the same integrated hard X-ray and gamma-ray emission for dark matter 
models with a bb, t^t~ and dominant final state, as a function of the WIMP mass. While a bb final state, 

common to numerous supersymmetric dark matter models, produces a soft electron-positron spectrum, leptonic final 
states give rise to a harder emission. In turn, this corresponds to a larger IC flux in the hard X-ray band, and to a 
suppressed gamma-ray flux. To compute the multi-wavelength emission from WIMP annihilation, we use the same 
setup as the one outlined in [13]. Notice that the ratio of the gamma-ray to hard X-ray luminosity from WIMP 
annihilation is largely independent of the particular setup chosen for the dark matter profile and for substructures 
(see SecHHI). 

We find that the generic expectation is for the ratio of gamma-ray to hard X-ray fiux to be larger for cosmic rays 
and suppressed for dark matter. The suppression can be as large as a factor 100, for final states producing a hard 
electron-positron spectrum, such as the lightest Kaluza-Klein particle of Universal Extra Dimensions (UED) [i^]. As 
we found in the spectral analysis and with the hardness ratio approach, the hardest cases to differentiate dark matter 
from cosmic rays are those with a very light dark matter particle. Comparing the left and the right panel of Fig. 1111 
we also find that smaller values of the magnetic field enhance the difference in the gamma-to-hard X-ray luminosity for 
cosmic rays and for dark matter, while larger values tend to give a more blurry picture. While the assessment of the 
average intra-cluster magnetic field (via e.g. Faraday rotation measurements) might only give an order of magnitude 
estimate, we find that for expected values for the magnetic field and for dark matter particle models with electro- weak 
scale particles, the gamma-ray to hard X-ray ratio technique proposed here would give a rather robust handle. Notice 
that in theoretically favored dark matter particle models, such as supersymmetry, the ratio shown in fig. [11] would 
typically lie between the black solid line (bb) and the green dashed line (t+t^). Other dark matter models, such as 
UED would have an even lower such ratio. 

In short, we showed that the ratio of the integrated gamma-ray to hard X-ray fiux in galaxy clusters can be used 
as a diagnostic for the discrimination of the origin of non-thermal phenomena, specifically astrophysical cosmic rays 
from dark matter annihilation. The generic expectation is that a dark matter induced signal would produce a brighter 
hard X-ray emission as opposed to cosmic rays, for a given detected gamma-ray fiux. 



IV. GAMMA RAY EMISSION FROM SELECTED GALAXY CLUSTERS AND GROUPS 



In this Section, we ask the question what are the best targets for the detection of gamma ray emission from clusters. 
We consider two large catalogs o f gal axy clusters and groups, the HIGFLUCS Catalog |llOl | and a subset of the catalog 



produced by the GEMS project [lllj . We describe below the assumptions we make to predict the gamma ray emission 
from cosmic rays and from dark matter for the objects considered in the two catalogs. In order to compare all objects, 
we make the same set of assumptions as far as both cosmic rays and dark matter are concerned: for instance, we 
assume the same ratio of gas to cosmic ray energy density {Xp = 0.01, i.e. lower by a factor 10 than what we 
considered before for Coma) for all clusters and groups. In modeling the dark matter halo, we assume the same 
fraction of mass in substructures versus the host halo (/s = 0.5) and the resulting = 7 x 10^ as for the "Boosted" 
setup (the fluxes corresponding to the "Smooth" case would have been roughly a factor 25 smaller, although this 
number depends on the cluster/group under consideration). 

To compute the cosmic ray emission, the first step is to extract the electron scaling density rie from the X-ray data. 
We assume that the intra-cluster medium (ICM) density can be described by a beta model: 

/ ^2 X -3/3/2 

pM-po 1 + ^ . (11) 



We follow here the approach outlined in [112l |. which assumes that the ICM is isothermal with temperature Tx, and 
we compute the scaling density of electrons Ue as 



2ntXeAiTx)RlFiRx) 



where 



which corresponds to the assumption of fully ionized plasma with 30% solar abundances [ll3j, Lx is the X-ray 



Pe,p = 4t ' ^ith He 1.167, ^1H = 1.400 (13) 

ne,H[r) rup 
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luminosity, Rx is the X-ray detection radius, and the function F is defined as [ll2| | 



F{R) 



ds 



{l + s' + {R/R,y) 



2n1-3/3 



2n1-3/3 



For the radiative coohng coefficient A(T) we assume the parametrization of jll4l |. i.e. 

A(r) = Ci(fcT)° + C2ikTf + C3, 

with 

Ci = 8.6 X 10"^^ erg cm^ s~^ keV"" 
C2 = 5.8 X 10^^" erg cm^ s"^ keV"'' 
C3 = 6.3 X IQ-^^ erg cm^ s"^ 
a = -1.7, /3 = 0.5. 



(14) 



(15) 



(16) 



We then use the cluster X-ray temperature Tx, the beta profile parameters Rc and (3, the cluster or group rcdshift 
and -R200 to compute the gamma-ray flux from cosmic rays. We assume an injection spectral index ap ~ 2.5 (this 
parameter is unimportant here, as it just re-normalizes all the cosmic-ray induced gamma-ray fluxes but not the 
relative emission), and we normalize the cosmic ray energy density to a fraction Xp = 0.01 of the gas energy density. 
For the cosmic ray source bias exponent /3p we consider two opposite cases: a "smooth" case, where ncn oc uh, hence 
Pp = 1 (see Sec. Ill A( ). and a "cuspy" case, where instead f3p = —0.5. We then integrate the gamma-ray emission 
above 0.1 GeV and out to i?2oo- Notice that this differs from the fluxes reported in Tab. Illli where we integrated the 
gamma-ray flux over the entire angular region corresponding to our Fermi simulations (10.5°). 

As far as the dark matter annihilation signal is concerned, we again use X-ray data to infer the dark matter density 
profile, assuming, as explained above, a Navarro- Frenk- White density distribution [t^I- We assume the substructure 
setup outlined in Sec. Ill Bl and in [sHi (to which we refer the reader for further details), and the structure formation 
model of [s^l- Here, we assume again that the radial density distribution of substructure follows the same profile as 
Eq. (O, but with a biased length scale a' ~ 7a, see e.g. [zEIzll; we also assume an average bias in the concentration 
of subhalos versus host halos at equal mass {cg) / (c^n-) = 4. Lastly, we assume that the fraction of mass in subhalos 
over the host halo mass fg = 0.5 (this setup corresponds to the Boosted or "B" case of sec. Ill Bp . We then determine 
the scale radius a and the scaling density po.DM in the following two ways. First, we take the (i?200i -^200) ^^nd 
(i?50o, pairs as determined from X-ray data, and solve for a and po,dm- Alternatively, we use the model of [s^l 

to relate the concentration of a given cluster to its mass and redshift (in practice for our low-redshift samples the 
latter does not matter significantly), and only make use of one of the (i?, M) pairs from X-ray data. Remarkably, we 
find that the two procedures yield very similar results in the final dark matter induced gamma-ray fiux, to the level 
of better than 10%. 

On the particle dark matter side, we assume a particle mass of 40 GcV, a pair annihilation cross section (av) = 
6 X 10~^^ cm^/s and 23 gamma rays per dark matter annihilation event above 0.1 GeV, as Monte Carlo simulations 
indicate is the case for a bb annihilation final state, ubiquitous e.g. in supersymmetric models j4l| . This particle 
DM model corresponds to the low mass (L) model of Sec. IIIBI Once the DM setup is specified, the gamma ray fiux 
from dark matter annihilations is then simply the integral over the line of sight of the number density squared of 
dark matter particles, times the annihilation rate. The resulting fiuxes we report, (f>o, can easily be rescaled for other 
masses, pair annihilation cross sections and annihilation final states: 



/40 GeVV 
V'TiwiMP/ V6 X 10"^^ cm^/s 



(av) 



' j^E>0.1 GcV ' 

23 



(17) 



In the case of the GEMS catalog |lll| . the group masses were not indicated. We need however that information 
to reconstruct the dark matter density profile. We thus assume that the ICM is in hydrostatic equilibrium and 
isothermal, and set |llOl | 



Af(< r) 



3kTx r^(3 f 1 



^mpGN + Rc , 

where nip is the proton mass, /i ~ 0.61 is the mean molecular weight and Gn is the gravitational constant. 



(18) 
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FIG. 12: Upper panel: predicted gamma-ray flux from dark matter annihilation in the extended HIGFLUCS catalog, with the 
assumptions for the dark matter profile specified in the text, and for a particle mass of 40 GeV, a bb dominant annihilation final 
state and a pair annihilation cross section of 6 x lO"'^^ cm^/s. Lower panel: X-ray data-driven predictions for the gamma-ray 
emission from cosmic rays in clusters and groups in the extended HIGFLUCS catalog. The red line refers to a cuspy spatial 
distribution for the cosmic ray sources (/3p = —0.5), while the black line to a smooth distribution (/?p = 1). 



A. Gamma-ray Emission from Clusters: the HIGFLUCS Catalog 



We applied the proc edure described above to the HIGFLUCS Catalog, including the clusters and groups from the 
Extended Sample [llOl |. The HIGFLUCS Catalog includes candidates from several input catalogs, and it includes 63 
clusters featuring an X-ray flux in the 0.1-2.4 keV range larger than 2.0 x 10~^^ ergs s~^ cm~^, with galactic latitude 
b > 20.0° and outside two excluded areas towards the Magellanic clouds and the Virgo cluster. In addition, we also 
include the Extended Sample, with 43 more clusters, bringing the total number of clusters to 106. Wc collect the flux 
ID numbers, names, predicted gamma-ray fluxes, ranking and ratio of cosmic ray to dark matter signal in Table IVll 
IVlIl and lVlIIl in the Appendix. 

We show the flux of gamma rays from dark matter annihilation in the upper panel of Fig. 1121 We flnd that the 
eight clusters with the expected largest dark matter induced gamma ray flux are the nearby large clusters Perseus, 
Coma, Ophiuchus and Abell 1060, 3526 and 3627. In addition the two groups M49 and Fornax also have very large 
fluxes, the latter in particular giving the largest one in the entire sample. The scatter in flux over the HIGFLUCS 
Catalog ranges over more than two orders of magnitude, but for our nominal choices for the dark matter particle 
properties and density distribution the typical flux from dark matter from low-redshift clusters lies typically between 
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FIG. 13: The ratio of the gamma-ray emission from dark matter and from cosmic rays for the clusters and groups in the 
extended HIGFLUCS catalog. The red and black lines refer to a cuspy and a smooth cosmic ray source distribution, as in 

Fig. \m 



10-1" and 10-^ g-i. 

The lower panel of Fig. [T^ shows our results for the cosmic-ray induced gamma-ray flux. We show the results 
for both a smooth (/3p = 1) and a cuspy (/3p = —0.5) primary cosmic ray source distribution. The cuspy source 
distribution typically boosts the flux by around one order of magnitude, and is expected e.g. in clusters with bright 
active galactic nuclei. The actual cosmic ray flux is expected to be somewhere in between the red and the black line. 
The clusters expected to be brightest in gamma-rays from cosmic ray interactions include Perseus, Coma, Ophiuchus, 
AbcU 3526, 3571, 3627 and 2319, as well as the bright cooling flow clusters 2A0335 and PKS0745. 

We remark that although our data-driven analytic approach is different and new with respect to other attempts at 
predicting which clusters are brighter in gamma rays, we substant ially agree with previous analyses, including psj 
which uses the ratio of the cluster mass over distance squared, and |l07l |. which makes use of the results of numerical 
simulations to assess scaling relations that then are used to predict the gamma-ray emission. We naturally agree with 
[3^, which uses a similar approach to ours. We note that these previous analyses predict only the gamma ray flux 
from cosmic rays and neglect dark matter annihilation. 

Fig. [T51 shows the ratio of the dark matter to cosmic ray gamma-ray flux for the extended HIGFLUCS sample. As 
in Fig.[T2]wc use the black line to indicate fluxes corresponding to a smooth cosmic ray source profile and a red line for 
the cuspy profile. In either approach, we notice several outliers, featuring a very high dark matter induced emission 
compared to the cosmic ray contribution. Candidates include NGC 5846, 5813, 499, 5044 and 4636, as well as Abell 
1060, 0548w, Fornax and M49. Interestingly, all these candidates are nearby poor clusters or galaxy groups. Even 
more interestingly, some of them (like M49 and Fornax) feature some of the largest predicted gamma-ray emissions 
from dark matter annihilation, making them particularly interesting candidates for Fermi observations. Even with 
extreme assumptions on the cosmic ray emission (/3p = —0.5), we find that the dark matter signal could be a factor 10 
larger than that of the cosmic ray induced flux. Searches for gamma-ray signals from these promising set of objects 
will be extremely valuable even with only one year of Fermi data. 

We show in Fig. [TJ] the correlation between the cluster mass to distance squared ratio and the predicted gamma-ray 
flux from cosmic rays and dark matter. We flnd a very tight correlation for the dark matter emission, which mainly 
depends on what we assumed for the reconstruction of the dark matter density proflles. On the other hand, the 
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FIG. 14: Correlation between the gamma-ray fluxes (from cosmic rays and darli matter) and the ratio of mass to distance 
squared, for objects in the HIGFLUCS catalog. 



scatter in the beta model parameters induce a significant scatter for the cosmic ray induced gamma-ray emission, 
with significant outhers both in excess and in deficit. 

Fig.[T5]shows a correlation betweeir the predicted gamma-ray luminosity from cosmic rays and the X-ray luminosity 
over distance squared for clusters in the extended HIGFLUCS sample. Although there are a few outliers, the corre- 
lation is rather tight, and close to linear. The dark matter induced emission, instead, does not show any significant 
correlation with the X-ray lumiirosity. Also, we do not find any other strong correlations between the predicted 
gamma-ray fluxes and X-ray related quantities. 

As a side comment, we remark that two of the three clusters for which |ll5l | tentatively associates an unidentified 
EGRET source with radio sources in the NRAO VLA aird Westerbork Northern sky survey catalogs appear also in 
our list, namely Abell 85 and 1914. While we predict Abell 85 be quite luminous in cosmic-rays-induced gamma-ray 
emission (ranking 12th to 16th out of 130, depending on the cosmic ray model), Abell 1914 is not predicted to be 
particularly bright. In addition, both clusters have a low dark matter to cosmic ray gamma-ray luminosity ratio, 
ranking respectively 124th and 87th. A dark matter interpretation for the tentative gamma-ray emissioir from these 
clusters seems therefore disfavored. 



B. Gamma-ray Emission from Groups: the GEMS Catalog 



Given the results of the previous Section, we decided to consider an alternate sample of nearby groups to (1) check 
whether our predictions with the HIGFLUCS sample depeird on the X-ray data analysis performed in [llO| | and (2) 
identify other potentially promising candid ates f or gamma ray searches for particle dark matter. We thus considered 
the catalog provided by the GEMS project [llll |. We excluded from our sample the following: 



1. groups with an X-ray flux less than 3cr above the background level (U-sample [lll| ) 



2. groups where the detectable extent of group emission is Text < 60 kpc, i.e. where the emission ap pear s to be 
associated with the halo of an individual galaxy instead of being geiruinely intergalactic (H-sample |lll| ) 



3. groups for which some of the relevant X-ray information was not available 
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FIG. 15: Correlation between the gamma-ray fluxes from cosmic rays and the ratio of X-ray luminosity to distance squared, 
for objects in the HIGFLUCS catalog. 



4. groups already included in the HIGFLUCS sample 

After the cuts described above, the GEMS sub-sample we used was limited to 24 groups. In passing, we notice that 
for the groups in the last item we obtained a remarkable agreement between the predictions for the gamma-ray fluxes 
we obtained with the HIGFLUCS and with the GEMS data. 

We show in Fig. [11] our results for the dark matter and cosmic ray induced gamma-ray fluxes (left) and the dark 
matter to cosmic ray ratio (right), with the same conventions as we used in Fig. [12] and [131 We find that the 
gamma-ray emission associated to the groups in the GEMS sample is generieally rather faint, with only one group 
(NGC1407) predicted to have a dark matter signal above 2 x 10"^'' cm"^ and only three groups (HCG62, NGC4073 
and NGC1407) that might have a gamma-ray emission from cosmic rays at the level of ~ 10~^° cm~^ s~^. More 
importantly, though, we confirm with the GEMS sample that nearby groups have a potentially very large ratio of dark 
matter to cosmic ray gamma-ray emission. We find in fact several candidates with a 100 times larger dark matter 
than cosmic ray gamma-ray emission, assuming a smooth cosmic ray source profile (see the right panel of Fig. I16p . 



V. DISCUSSION, SUMMARY AND CONCLUSIONS 



The detection of gamma rays from clusters of galaxies might be a milestone in the scientific legacy of the Fermi 
Gamma-ray Space Telescope. Clusters are known to host a variety of high-energy phenomena that could fuel cosmic 
rays, producing, in turn, gamma rays as a result of collisions with the ICM gas. Being the largest bound dark matter 
structures, it is also reasonable to expect that clusters feature a significant gamma-ray emission from the annihilation 
of dark matter particles. The scope of the present theoretical study was to investigate how to tell apart these two 
potential mechanisms of gamma-ray production with data from Fermi. 

One potential source of gamma-rays from cluster regions that we have not considered here are the bright, central 
AGN known to be present in some clusters (NGC 1275 in Perseus and M87 in Virgo), which could inhibit our ability 
to detect faint extended gamma-ray emission. However, this source of contamination is not expected from some of the 
best candidate, nearby clusters, like the Coma cluster, where strong radio galaxies are not observed. The best method 
to detect diffuse gamma-ray emission, from cosmic rays or dark matter annihilation, will be to concentrate on those 
clusters lacking bright AGN. For clusters with central AGN, the expected extended nature of the cluster gamma-ray 
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FIG. 16: Gamma-ray fluxes from dark matter and cosmic rays (left) and their ratio (right) for our sub-sample of groups in the 
GEMS catalog. 



emission may reveal this component, and the well-known multiwavelength (radio and X-ray) AGN spectra will allow 
modeling of the contribution of these sources to the gamma-ray emission. 

By making simple analytical assumptions on the cosmic ray spectra and source distribution and on the dark matter 
particle properties and density distribution, wc proposed a set of various benchmark models for both gamma-ray 
production mechanisms under investigation. We believe the set of models we considered is representative of the 
variety of possibilities one can realistically expect to encounter in clusters of galaxies. We then proceeded to simulate 
the expected 1 and 5 years gamma-ray signals with the Fermi Science Tools, and analyzed our results. 

We summarize below the main results of the present theoretical study. 

• We find that for cosmic rays the absolute signal to noise peaks at around one degree, while for dark matter 
it can peak at larger regions of interest, possibly as large as 3 degrees. Since the Fermi instrument response 
function forces us to consider relatively large angular regions to include low energy photons, we conclude that 
the best region of interest for gamma-ray studies of cluster of galaxies is around 3 degrees. 

• The spectral analysis of the simulated signal for our benchmark models shows that gamma ray emission from dark 
matter annihilation with a relatively large dark matter particle mass (ttiwimp > 50 GeV) can be distinguished 
from a cosmic ray spectrum even for fairly faint sources. Distinguishing a cosmic ray spectrum from a low dark 
matter particle mass appears to be more challenging, and would require deep data and/or a bright source. 

• We defined optimal energy bands for a simple hardness ratio, (H - S)/(H -I- S), estimate of the spectrum. The 
energy bands we propose are S = 0.15 < E < 0.7 GeV and H = E > 0.7 GeV. The hardness ratio is correlated 
to the nature of the emission and is predicted to be negative for a cosmic ray emission and positive for a dark 
matter annihilation signal. Similar to the full spectral analysis, the hardness ratios for a simulated low mass 
dark matter particle model and a hard cosmic ray model are similar within the statisitcal errors. However, it is 
possible to use this simple two energy band ratio to distinguish cosmic ray and high mass dark matter particle 
models. 

• The level of the gamma-ray emission from clusters produced by cosmic rays can be comparable to that from 
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dark matter annihilation. In the case of a mix of the two emissions, we find that for bright enough sources, if the 
dark matter contribution to the cluster flux is significant and the particle mass is not very low {mwiMP > 40 
GeV) the presence of a dark matter component can be seen even in the presence of a significant gamma-ray flux 
from cosmic rays. However, tight constraints on the model parameters {mwiMP and ap) may be problematic 
unless the dark matter particle mass is known. 

• Our cluster gamma-ray simulated emissions appear, after data reduction, as extended rather than point sources, 
with extensions which depend on the spatial models and on the emission mechanism. However, determining the 
spatial distribution with Fermi will be a challenging task requiring an optimal control of the backgrounds. 

• We showed that the ratio of the integrated gamma-ray to hard X-ray flux in galaxy clusters can be used as a 
diagnostics for the discrimination of the origin of non-thermal phenomena. The generic expectation is that a 
dark matter induced signal would produce a brighter hard X-ray emission as opposed to cosmic rays, for a given 
gamma-ray flux. 

• We presented X-ray data-driven predictions for the gamma-ray flux from 130 clusters and groups of galaxies 
in the HIGFLUCS and GEMS catalogs. We found that the clusters with the brightest gamma-ray emission 
from cosmic rays include the Perseus, Coma, Ophiuchus, Abell 3627 and Abell 3526 clusters; the most luminous 
clusters in dark matter emission are predicted to be the Fornax group, Ophiuchus, Coma. Abell 3526 and Abell 
3627. 

• We discovered that the objects with the largest dark matter to cosmic ray gamma-ray luminosity in our sample 
are groups and poor clusters. In particular, the highest ratios are associated to the groups NGC 5846, 5813 
and 499, M49 and HCG 22. Of these, M49 ranks overall 6th/130 in terms of predicted dark matter induced 
gamma-ray emission, NGC 5846 ranks 12th and NGC 5813 15th; Fornax also has a relatively high dark matter 
to cosmic ray gamma-ray flux, and is the brightest object in dark matter emission. All these objects are very 
promising candidates for a search for gamma-ray emission with the Fermi telescope that could potentially be 
related to particle dark matter: specifically, as shown in the Appendix, we predict, for the DM_HS setup, an 
integrated flux of photons above 0.1 GeV of 4 x 10~^ cm~^s~^ for Fornax and of 2 x 10"^ cm^^s^^ for M49, 
both around the projected Fermi point-source sensitivity for one year of data. With the same setup, we predict 
a flux of 0.8 X 10^^ cm~^s~^ and 0.5 x 10^^ cm^^s^^ for NGC 5846 and 5813, respectively, which would make 
them detectable or marginally detectable sources for Fermi within the anticipated mission lifetime. 
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Appendix 



We report here the X-ray data driven gamma-ray flux predictions for the HIGFLUCS catalog ftab. lVII and lVIIl plus 
members of the extended catalog fVlIip and for the GEMS catalog (tab. lIXp . The first column indicates the object ID 
used in our plots, the second one specifies the object name. The third column shows the predicted gamma-ray fiux 
from dark matter annihilation for the model setup DM_LB of sec. Ill Bl (featuring a particle mass of 40 GeV, a pair 
annihilation cross section of 6 x 10"'^^ cm'^/s, and 23 photons with an energy above 0.1 GeV per annihilation), and 
the fourth column the overall ranking (considering at the same time the extended HIGFLUCS catalog and the GEMS 
catalog). The next columns indicate the predictions for a cosmic ray origin, under the assumption of a smooth and 
of a cuspy source distribution, as well well as the overal ranking according to the corresponding emission. Lastly, the 
rightmost columns indicate the ratio of the dark matter to cosmic ray gamma-ray emission, with ranking. 

The fiux predictions quoted here can be easily rescaled for different particle DM models as: 



where (j)Q is the integrated gamma-ray fiux above 0.1 GeV we quote in the following tables and 0' is the corresponding 
re-scaled fiux. Similarly, for the cosmic-ray predictions, the scaling with the fraction of energy density in cosmic rays 
Xp (which we here set at 1%) is simply 



In the tables, we mark in red the top-ten objects in every category we consider (gamma-ray emission from DM 
annihilation, emission from cosmic-ray inelastic interactions, and ratio of the DM-to-cosmic-rays projected emission). 

In the figures, we also indicate with a horizontal dashed line an estimate of the high-latitude point source sensitivity 
of Fermi-LAT over 5 years (~ 1.3 x 10~^ photons per cm^ per s). We predict, with the assumptions we made here 
for the DM setup, that on the order of 10 clusters will produce a large enough DM annihilation emission of gamma 
rays. Over the anticipated Fermi lifetime of 10 years, 5 or so more clusters might also be detectable according to the 
present predictions. 




(19) 




(20) 
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ID Name DMXB DMXB, CR_S CR_S, CR_C CR_C, DMXB/CR_S DM_LB/CR_S, 









rank 




rank 




rank 




rank 


1 


A 0085 


213.1 


47 


110.3 


12 


1285.0 


16 


1.9 


124 


2 


A0119 


336.4 


30 


73.3 


25 


382.4 


47 


4.6 


86 


3 


A0133 


83.1 


98 


19.6 


74 


277.3 


61 


4.2 


92 


4 


NGC507 


141.3 


74 


11.8 


96 


85.5 


104 


11.9 


40 


5 


A0262 


329.2 


32 


74.7 


24 


446.4 


36 


4.4 


89 


6 


A0400 


216.6 


46 


24.5 


68 


132.6 


91 


8.8 


54 


7 


A0399 


199.5 


50 


56.6 


28 


418.2 


41 


3.5 


101 


8 


A0401 


182.2 


54 


87.9 


19 


801.0 


24 


2.1 


122 


9 


A3112 


89.8 


94 


28.0 


57 


561.3 


31 


3.2 


105 


10 


FORNAX 


4016.0 


1 


37.2 


43 


419.0 


40 


108.1 


16 


11 


2A0335 


174.4 


62 


48.5 


33 


1286.0 


15 


3.6 


99 


12 


IIIZw54 


245.3 


40 


12.0 


95 


135.8 


88 


20.4 


25 


13 


A3158 


198.1 


51 


54.5 


31 


496.8 


34 


3.6 


98 


14 


A0478 


136.5 


76 


56.4 


29 


1259.0 


17 


2.4 


116 


15 


NGC1550 


432.4 


22 


22.8 


71 


280.4 


60 


19.0 


27 


16 


RXO0422 


184.5 


53 


18.3 


78 


317.6 


55 


10.1 


50 


17 


A 3266 


410.8 


23 


103.9 


14 


858.9 


21 


4.0 


95 


18 


A0496 


244.5 


41 


90.4 


17 


1024.0 


20 


2.7 


112 


19 


A3376 


346.3 


29 


29.6 


53 


212.0 


71 


11.7 


42 


20 


A3391 


181.7 


56 


36.3 


45 


244.2 


69 


5.0 


83 


21 


A3395s 


375.0 


25 


25.7 


63 


237.3 


70 


14.6 


33 


22 


A0576 


373.6 


26 


31.4 


49 


293.0 


56 


11.9 


41 


23 


A0754 


572.4 


14 


46.6 


36 


774.2 


27 


12.3 


38 


24 


HYDR A- A 


125.1 


81 


37.2 


42 


782.2 


25 


3.4 


104 


25 


A1060 


1958.0 


7 


80.8 


23 


1137.0 


19 


24.2 


23 


26 


A1367 


732.9 


13 


81.2 


22 


467.1 


35 


9.0 


52 


27 


MKW4 


150.4 


70 


14.2 


89 


122.5 


94 


10.6 


45 


28 


ZwC11215 


158.5 


67 


27.6 


58 


275.9 


62 


5.7 


72 


29 


NGC4636 


1530.0 


9 


11.0 


99 


197.7 


72 


138.8 


14 


30 


A3526 


2137.0 


4 


273.0 


5 


3028.0 


5 


7.8 


58 



TABLE VI: Gamma-ray emission predictions for the groups and clusters in the HIGFLUCS catalog (first part). Fluxes are in 
units of 10"^^ph/(cm^ s) 
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ID Name DMXB DMXB, CR_S CR_S, CR_C CR_C, DMXB/CR_S DMXB/CR_S, 









rank 




rank 




rank 




rank 


31 


A1644 


182.1 


55 


68.1 


26 


359.5 


51 


2.7 


113 


32 


A1650 


134.7 


77 


31.5 


48 


367.7 


50 


4.3 


91 


33 


A1651 


99.3 


88 


31.3 


50 


400.9 


43 


3.2 


107 


34 


COMA 


2170.0 


3 


610.8 


2 


5073.0 


3 


3.6 


100 


35 


NGC5044 


481.1 


19 


17.0 


81 


375.1 


49 


28.4 


22 


36 


A 1736 


107.8 


85 


47.3 


34 


171.0 


81 


2.3 


118 


37 


A3558 


228.1 


43 


106.6 


13 


750.2 


28 


2.1 


121 


38 


A3562 


142.5 


73 


46.8 


35 


287.5 


57 


3.0 


109 


39 


A3571 


512.0 


16 


171.0 


9 


1929.0 


7 


3.0 


110 


40 


A 1795 


247.9 


38 


66.8 


27 


1546.0 


11 


3.7 


97 


41 


A3581 


199.9 


49 


19.2 


76 


284.4 


58 


10.4 


47 


42 


MKW8 


275.6 


34 


27.2 


59 


182.4 


75 


10.1 


48 


43 


A2029 


196.4 


52 


86.9 


20 


1787.0 


10 


2.3 


119 


44 


A2052 


154.9 


69 


38.3 


40 


548.6 


32 


4.0 


94 


45 


MKW3S 


145.2 


72 


26.2 


60 


438.6 


37 


5.5 


75 


46 


A 2065 


274.9 


35 


25.6 


64 


333.0 


53 


10.7 


44 


47 


A2063 


219.3 


44 


42.3 


38 


398.1 


45 


5.2 


80 


48 


A2142 


160.8 


66 


101.7 


16 


1333.0 


13 


1.6 


125 


49 


A2147 


237.8 


42 


103.8 


15 


375.9 


48 


2.3 


117 


50 


A2163 


83.1 


99 


40.3 


39 


508.0 


33 


2.1 


123 


51 


A2199 


445.8 


20 


88.2 


18 


1298.0 


14 


5.1 


81 


52 


A2204 


37.6 


120 


25.8 


62 


674.0 


29 


1.5 


128 


53 


A2244 


90.7 


92 


25.9 


61 


400.6 


44 


3.5 


102 


54 


A2256 


368.1 


27 


83.2 


21 


841.3 


22 


4.4 


88 


55 


A2255 


178.9 


58 


35.6 


46 


253.6 


66 


5.0 


82 


56 


A3667 


217.6 


45 


143.2 


11 


780.4 


26 


1.5 


126 


57 


SllOl 


74.0 


101 


11.5 


97 


339.7 


52 


6.4 


65 


58 


A2589 


175.6 


61 


25.2 


66 


283.4 


59 


7.0 


62 


59 


A2597 


60.7 


109 


12.9 


92 


427.5 


38 


4.7 


85 


60 


A2634 


331.5 


31 


34.6 


47 


177.0 


79 


9.6 


51 


61 


A2657 


167.9 


63 


28.6 


56 


246.7 


68 


5.9 


71 


62 


A4038 


257.0 


37 


49.0 


32 


592.2 


30 


5.2 


78 


63 


A4059 


180.9 


57 


30.7 


51 


427.3 


39 


5.9 


70 



TABLE VII: Gamma-ray emission predictions for the groups and clusters in the HIGFLUCS catalog (second part). Fluxes are 
in units of 10~^'^ph/(cm'^ s) 



34 



ID Name DMXB DMXB, CR_S CR^, CR_C CR_C, DMXB/CR_S DMXB/CR_S, 









rank 




rank 




rank 




rank 


64 


A2734 


90.2 


93 


17.0 


80 


134.3 


89 


5.3 


77 


65 


A2877 


440.0 


21 


23.3 


70 


146.9 


85 


18.9 


28 


66 


NGC499 


161.2 


65 


0.4 


123 


26.9 


110 


417.1 


4 


67 


AWM7 


1236.0 


10 


149.7 


10 


1819.0 


8 


8.3 


57 


68 


PERSEUS 


1940.0 


8 


1353.0 


1 


20230.0 


1 


1.4 


130 


69 


S405 


109.4 


84 


28.8 


54 


137.4 


87 


3.8 


96 


70 


30129^ 


1126.0 


11 


201.2 




1189.0 


18 


5.6 


73 


71 


A0539 


272.3 


36 


36.4 


44 


261.6 


63 


7.5 


59 


72 


S540 


138.4 


75 


12.1 


94 


133.7 


90 


11.4 


43 


73 


A0548w 


35.7 


122 


1.5 


111 


8.8 


117 


23.9 


24 


74 


A0548e 


99.9 


87 


16.4 


83 


85.7 


103 


6.1 


67 


75 


A3395n 


377.4 


24 


21.8 


73 


180.6 


78 


17.3 


30 


76 


UGC03957 


210.4 


48 


10.7 


100 


190.9 


73 


19.6 


26 


77 


PKS0745t 


83.0 


100 


55.9 


30 


1533.0 


12 


1.5 


127 


78 


A0644 


245.4 


39 


44.3 


37 


804.5 


23 


5.5 


74 


79 


S636 


510.1 


17 


37.9 


41 


165.8 


82 


13.5 


35 


80 


A1413 


50.5 


113 


15.8 


85 


249.6 


67 


3.2 


106 


81 


M49 


2009.0 


6 


2.8 


106 


149.3 


84 


710.0 


3 


82 


A3528n 


97.2 


91 


13.2 


91 


113.9 


98 


7.3 


60 


83 


A3528s 


54.2 


111 


22.4 


72 


115.8 


96 


2.4 


115 


84 


A3530 


157.1 


68 


12.7 


93 


89.4 


102 


12.4 


37 


85 


A3532 


163.4 


64 


25.2 


65 


189.0 


74 


6.5 


64 


86 


A1689 


46.8 


115 


15.9 


84 


387.7 


46 


2.9 


111 


87 


A3560 


89.1 


95 


18.6 


77 


90.6 


101 


4.8 


84 


88 


A 1775 


63.4 


108 


14.4 


88 


114.8 


97 


4.4 


90 


89 


A1800 


87.4 


97 


14.4 


87 


109.8 


99 


6.1 


69 


90 


A1914 


74.0 


102 


16.5 


82 


407.3 


42 


4.5 


87 


91 


NGC5813 


555.6 


15 


0.8 


118 


54.7 


107 


724.7 


2 


92 


NGC5846 


853.0 


12 


0.9 


115 


76.6 


105 


943.5 


1 


93 


A2151w 


107.6 


86 


10.7 


101 


126.2 


93 


10.1 


49 


94 


A3627''^ 


2081.0 


5 


598.1 


3 


3128.0 


4 


3.5 


103 


95 


TRIANGUL 


504.2 


18 


233.0 


7 


2018.0 


6 


2.2 


120 


96 


OPHIUCHU^ 


2497.0 


2 


459.8 


4 


8642.0 


2 


5.4 


76 


97 


ZwC11742 


118.5 


82 


19.5 


75 


261.1 


64 


6.1 


68 


98 


A2319 


346.5 


28 


240.5 


6 


1790.0 


9 


1.4 


129 


99 


A3695 


72.6 


104 


29.6 


52 


165.0 


83 


2.5 


114 


100 


TT7w1 OS 


125.4 


80 


24.0 


69 


127.8 


92 


5.2 


79 


101 


A3822 


87.8 


96 


28.7 


55 


171.1 


80 


3.1 


108 


102 


A3827 


176.9 


59 


24.8 


67 


320.3 


54 


7.1 


61 


103 


A3888 


99.2 


89 


11.5 


98 


254.7 


65 


8.6 


55 


104 


A3921 


97.4 


90 


15.7 


86 


180.8 


77 


6.2 


66 


105 


HCG94 


126.3 


79 


18.2 


79 


144.0 


86 


6.9 


63 


106 


RXJ2344 


112.2 


83 


13.3 


90 


181.2 


76 


8.4 


56 



TABLE VIII: Gamma-ray emission predictions for the groups and clusters in the Extended Sample of the HIGFLUCS catalog. 
Fluxes are in units of 10^^^ph/(cm^ s). Clusters at low galatic latitudes, |6| < 10 deg, are marked with f . These clusters will 
have higher gamma-ray backgrounds due to galactic emission. 
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ID Name DMXB DM_LB, CR^ CR_S, CR_C CR_C, DM_LB/CR_S DMXB/CR_S, 









rank 




rank 




rank 




rank 


1 


NOC383 


57.1 


110 


6.5 


104 


29.7 


109 


8.8 


53 


2 


NGC533 


36.9 


121 


2.0 


109 


17.2 


111 


18.2 


29 


3 


NGC720 


150.0 


71 


0.6 


121 


10.7 


116 


267.9 


7 


4 


NGC741 


42.9 


117 


1.1 


113 


12.1 


113 


39.4 


20 


5 


HCG22 


17.7 


126 


0.05 


130 


0.5 


130 


383.6 


5 


6 


NGC1407 


321.0 


33 


2.3 


108 


93.3 


100 


141.5 


13 


7 


NGC1587 


65.1 


107 


0.3 


124 


3.9 


127 


190.6 


10 


8 


NGC2563 


33.7 


123 


2.3 


107 


11.5 


115 


14.6 


34 


9 


HCG42 


40.0 


118 


0.2 


126 


7.3 


119 


206.3 


g 


10 


NGC3557 


18.0 


125 


0.1 


128 


5.1 


123 


132.9 


15 


11 


NGC3607 


69.1 


105 


1.0 


114 


5.9 


122 


66.0 


18 


12 


NGC3665 


49.8 


114 


0.1 


127 


2.7 


128 


352.3 


6 


13 


NGC4065 


19.6 


124 


1.6 


no 


7.1 


120 


12.4 


36 


14 


NGC4073 


39.7 


119 


9.5 


103 


66.5 


106 


4.2 


93 


15 


NGC4261 


176.0 


60 


10.2 


102 


46.6 


108 


17.2 


31 


16 


NGG4325 


14.3 


128 


1.4 


112 


14.6 


112 


10.6 


46 


17 


NGC4589 


129.4 


78 


0.9 


116 


11.6 


114 


148.6 


12 


1 8 


HGGfi9 


fiS 1 


-LUU 


O. 1 


1 (Ti 


117.2 


yu 


12.0 


oy 


19 


NGC5129 


12.8 


129 


0.9 


117 


7.3 


118 


15.0 


32 


20 


HCG67 


10.4 


130 


0.1 


129 


1.7 


129 


170.6 


11 


21 


HCG68 


54.2 


112 


0.6 


120 


4.8 


125 


94.6 


17 


22 


HCG90 


44.9 


116 


0.7 


119 


6.0 


121 


60.2 


19 


23 


IC1459 


73.8 


103 


0.3 


125 


4.6 


126 


239.2 


8 


24 


HCG97 


17.6 


127 


0.5 


122 


4.8 


124 


36.6 


21 



TABLE IX: Gamma-ray emission predictions for the groups in the GEMS catalog. Fluxes are in units of 10"^^ph/(cm^ s) 



